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Q , Abstract. In previous work we introduced and studied a function i?(a+, a_, c; v, v) that 

is a generalization of the hypergeometric function 2-P1 and the Askey-Wilson polynomials. 
When the coupling vector c 6 C 4 is specialized to (6,0,0,0), b e C, we obtain a func- 
tion Tl(a+, a_, b; v, 2v) that generalizes the conical function specialization of 2F1 and the 
q-Gcgcnbaucr polynomials. The function 1Z is the joint eigenfunction of four analytic diffc- 
rence operators associated with the relativistic Calogero-Moser system of A\ type, whereas 
\ the function R corresponds to BCi, and is the joint eigenfunction of four hyperbolic Askey- 

Wilson type difference operators. We show that the 7^-function admits five novel integral 
representations that involve only four hyperbolic gamma functions and plane waves. Taking 
their nonrelativistic limit, we arrive at four representations of the conical function. We also 
show that a limit procedure leads to two commuting relativistic Toda Hamiltonians and 
two commuting dual Toda Hamiltonians, and that a similarity transform of the function 1Z 
converges to a joint eigenfunction of the latter four difference operators. 
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1 Introduction 

This article may be viewed as a continuation of our previous work on a 'relativistic' generaliza- 
tion R of the Gauss hyper geometric function 2-F1, introduced in [JJ. The latter paper and two 
later parts in a series [2j [3] will be referred to as I, II and III in the sequel. The definition of 
the .R-function in I is in terms of a contour integral, whose integrand involves eight hyperbolic 
gamma functions. (We review this in Section O cf. (j2.ip - (j2.5p .) 

In recent years, van de Bult [1] tied in the R- function with the notion of modular double of 
the quantum group U q (sl(2, C)), as defined by Faddeev (5j. As a spin-off, he obtained a new 
representation of the i?-function. Also, van de Bult, Rains and Stokman [6] have shown (among 
other things) that the 8-variable i?-function 

R(a + ,a-,c;v,v), a+, a_,v,v £ C, a + /a_ ^ (—00, 0], c £ C 4 , (1.1) 

can be obtained as a limit of Spiridonov's 9- variable hyperbolic hyper geometric function [7]. 
Their novel viewpoint leads to a third representation for the i?-function. (See Proposition 4.20 
and Theorem 4.21 in (6j for the latter two representations.) 

In this paper we are concerned with a 5-variable specialization of the i?-function, defined by 

K(a+,a-,b;x,y) = R{a + ,a., (6, 0,0,0); x,y/2). (1.2) 

Suitable discretizations of this function give rise to the g-Gegenbauer polynomials, whereas 
discretizations of the i?-function yield the Askey-Wilson polynomials, cf. I; moreover, the non- 
relativistic limit of the 1Z- function yields the conical function specialization of 2-F1. Hence it 
may be viewed as corresponding to the Lie algebra A±, whereas the R- function can be tied in 
with BCi. 

The key new result of this paper concerning 1Z consists of the integral representation 

f^ G(2ib - ia) f G(z + (x- y)/2 - ib/2)G(z - (x - y)/2 - ib/2) 
1 ' X ' V) V 2vr G{ib - iaf U Z G(z + (x + y)/2 + ib/2)G{z - (x + y)/2 + ib/2) ' 1 ' ' 

Here and throughout the paper we use parameters 

a = 27r/a + a_, a = (a + + a_)/2, (1.4) 

G(a + , a_; z) is the hyperbolic gamma function (cf. Appendix|A]), and the dependence on a + , a_ 
is suppressed. (We shall often do this when no confusion can arise.) Furthermore, in fll .3 j) we 
choose at first 

(a + , a_, b, x, y) G (0,oo) 2 x (0,2a) x M 2 . (1.5) 
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By contrast to the previous three integral representations following from I, [5] and [6], the 
integrand in (|1.3|) involves only four hyperbolic gamma functions. We also obtain several closely 
related representations that involve in addition plane waves, cf. (|3.47p - (|3.5ip . As will transpire in 
Section [3l upon using the first one (jl.3p of these novel representations (which we dub 'minimal' 
representations) to introduce the 7£-function, it is possible to rederive in a more transparent 
and self-contained way a great many features that also follow upon specialization of the R- 
function theory, developed not only in I, II and III, but also in our later papers [8] and [9]. 
Moreover, special cases and limits of the 7£-function are far more easily obtained from the 
minimal representations than from the original integral representation of I or from the alternative 
representations following from [I] and [6]. (The integrands of these earlier representations involve 
at least eight hyperbolic gamma function factors.) 

A survey of the results of I— III and [8] can be found in [10], but the definition (jl.2D of the 
Ai-analog of the (BC\) R- function dates back to the more recent paper [9]. In Section [21 we 
review in particular the pertinent results from [9]. However, we have occasion to add a lot more 
information that follows by specializing previous findings concerning the i?-function and related 
functions to their A± counterparts. This includes the asymptotic behavior and Hilbert space 
properties obtained in II and III, resp., which are adapted to the A± setting in Subsection 12.21 
and the parameter shifts obtained in [8], which we focus on in Subsection [231 Moreover, in (I2.27j) 
we detail the connection of the renormalized function 



to the Ai type functions M(ma+ + ?ia_; x,y), m,n E Z, which featured in our previous pa- 
pers [TT] and [12]. We present the proof of (|2.27p in Subsection 12.31 together with various 
corollaries. 

Altogether, Section [2] invokes a considerable amount of information from our previous work. 
We have attempted to sketch this in such a way that the reader need only consult the pertinent 
papers for quite technical aspects (in case of doubt and/or inclination, of course). Even so, it 
is probably advisable to skim through Section [21 at first reading, referring back to it when the 
need arises. 

By contrast, Section [3] (combined with Appendices A and C) is largely self-contained. Its 
starting point is a hyperbolic functional identity that first arose as a specialization of elliptic 
functional identities expressing the relation of certain Hilbert-Schmidt integral kernels to the 
elliptic BC\ relativistic Calogero-Moser difference operators introduced by van Diejen [13] . We 
need not invoke these identities (which can be found in [14] . cf. also [E])> since the relevant 
hyperbolic version is quite easily proved directly. The key point is that the hyperbolic identities 
can be rewritten in terms of two pairs of hyperbolic ^4-i-type relativistic Calogero-Moser diffe- 
rence operators A±(bj-,x), j = 1,2, with distinct couplings b±, &2- The difference operators are 
given by 



1l r (a + ,a_,b;x,y) 



G(ib — ia) 
G(2ib - ia) 



Tl(a + ,a-,b;x,y) 



(1.6) 



A s (b;x) 



s$(x — ib) 



s s (x) 




sg(x + ib) 
s s (x) 



5 = + 



(1.7) 



Here, the translations are defined on analytic functions by 



(1.8) 



Also, throughout this paper we use the abbreviations 



sg(z) = s'mh.(ir z / as) , c$(z) = cosh^z/a^), es{z) = exp(irz/as), 5 = +, — . (1.9) 
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Choosing b\ = b, &2 = 0, an auxiliary function B(b;x,y) can be defined that satisfies the 
eigenvalue equations (at first under certain restrictions on the £>-arguments) 

A s (b;x)B(b;x,y) = 2c 5 (y)B(b;x,y), 5 = +,-. (1.10) 

More specifically, the function B(b;x,y) is the Fourier transform of the hyperbolic kernel func- 
tion, which is a product of four hyperbolic gamma functions. 

When we write the integrand of the integral defining B as a product of two factors that involve 
only two hyperbolic gamma functions, we can use the Plancherel relation and the explicit Fourier 
transform formula for factors of this type (derived in Appendix [C]) to obtain two new integral 
representations for B. In particular, this leads to a function C(b;x,y) given by 

r(b-rv)= [° f d~ G{z + {x ~ y)/2 ~ lb,2)G{z - {x ~ y)/2 ~ lh/2) am 

{ ' ' U1 ~ V 2tt J r G(z + (x + y)/2 + ib/2)G(z - (x + y)/2 + ib/2) ' 1 ' } 

Comparing (jl.lip . (|1.3|) and (|1.6p . we read off 

K r (b; x, y) = G(ib - ia)C(b; x, y). (1.12) 

However, we need a further study of the C-function (jl.lip to arrive at a proof of this relation 
to the function lZ r , as defined originally by fjl .6|) and (jl.2p . Indeed, as already alluded to 
below (|1.5|) . we can use (jl.lip as a starting point to derive many features that C and 7Z r have 
in common. 

In particular, the general analysis in Appendix B of I can be applied to the integral on 
the r.h.s. of (jl.lip . which yields a complete elucidation of the behavior of C(b;x,y) under 
meromorphic continuation. Moreover, via the AAEs (analytic difference equations) (jl.lOp and 
the manifest invariance of C under interchanging x and y, it follows that C(b;x,y) is a joint 
eigenfunction of the four AAOs (analytic difference operators) 

A + (b;x), A_(6;x), A+(b;y), A_(&;y), (1.13) 

with eigenvalues 

2c+(y), 2c_(y), 2c+(x), 2c_(x), (1.14) 

resp. This is also the case for lZ r (b;x,y) and, moreover, the equality (11 . 12f) can be shown for 
the special case y = ib by a further application of Appendix [Cj The general case then follows 
by a uniqueness argument already used in Subsection 12.31 
We reconsider the special 5-values 

6 mn = ma + + na„, m, n£Z, (1-15) 

in Subsection 14. lj inasmuch as they satisfy b mn G (0,2a). Indeed, the new Fourier transform 
representations in Section [3] are only well defined for b G (0,2a), but they can be explicitly 
evaluated by a residue calculation when b is of this form. The key point is that the G-ratios 
in the integrand can then be written in terms of the hyperbolic cosines c±(z) by using the G- 
AAEs (jA.2p . In principle, this yields again the functions M(b mn ;x,y) from [11], but we have 
not tried to push through a direct equality proof (as opposed to appealing to uniqueness). 

Subsection 14.21 deals with the nonrelativistic limit. Specializing the results of I yields the 
hyper geometric function in terms of which the conical function can be expressed (cf. Chapter 14 
in [16]). The five minimal representations (j3.47p - (j3.5ip of the 1Z- function lead to four repre- 
sentations (j4.45p ~ (j4.48p of the limit function. Rewriting them in terms of the conical function, 
three of these can be found in the literature (by looking rather hard). This is reassuring, since 
just as in I we were not able to get rigorous control on the nonrelativistic limits. 
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In order to describe the results of Section [5j we begin by recalling that in our paper [T7] we 
arrived at relativistic nonperiodic Toda iV-particle systems by taking a limit of the relativistic 
hyperbolic Calogero-Moser iV-particle systems. In this limit the self-duality of the latter is 
not preserved, inasmuch as the dual commuting Hamiltonians have a very different character 
from the defining Hamiltonian and its commuting family. Specialized to the present context, 
this limit can be used to obtain a joint eigenfunction of two Toda Hamiltonians H±(rj;x) and 
two dual Toda Hamiltonians H±(r]; y), with the real parameter 77 playing the role of a coupling 
constant. 

The limit transition proceeds in two stages. The first step is to set 

b = a -i~f, 7GM. (1.16) 

At the classical level the analogous 6-choice still yields real-valued Hamiltonians with a well- 
defined self-dual action-angle map and scattering theory |18j . Correspondingly, the four reduced 
N = 2 quantum Hamiltonians at issue here are still formally self-adjoint for this 6-choice. (They 
are similarity transforms of the AAOs (11. 13ft with a weight function factor.) Moreover, restricting 
attention to 

(a + ,a-,x,y) £ (0,oo) 2 x M 2 , (1.17) 

their joint eigenfunction remains real- valued, although this reality property is no longer manifest: 
It hinges on a symmetry property under taking b to 2a — b, which translates into evenness in 
the parameter 7. 

The next step is to substitute 

x^x + A, 7-^77 + A, (1-18) 

and take A to 00. In this limit the Hamiltonians and their joint eigenfunction converge, whereas 
the dual Hamiltonians must be multiplied by a factor e<y(— A) to obtain a finite limit. This can 
be understood from their A-dependent eigenvalues 2c$(x + A) following from the x-shift (|1.18p . 
cf. (|1.14p . Indeed, after multiplication by this renormalizing factor the eigenvalues have the 
finite limits eg(x), 5 = +,—. 

The five representations of the 7£-function give rise to four representations of the relativistic 
Toda eigenfunction F T (rj;x,y), namely f|5.25j) . (|5.26p . (|5.32p and (|5.33p . Suitably paired off, 
however, these different formulas express real-valuedness with (|1.17p in effect. Taking this 
into account, we wind up with two essentially different representations that are intertwined 
via the Plancherel formula for the Fourier transform. The key formula involved here is derived 
in Corollary IC. 21 

The results just delineated can be found in Subsection 15.11 In Subsection 15.21 we first study 
the asymptotic behavior of F T (r]; x, y) for x — > ±00 and y — > 00. We then clarify the analyticity 
properties of F T (rj;x,y) by introducing a similarity transform %{x — rj,y). Using the four 
representations (|5.54p - (l5.57p of the latter, we show that the function T~L(x,y) is holomorphic for 

(x,y)eC 2 . 

Subsection 15.31 deals with the joint eigenfunction properties of F T (rj;x,y) and its similarity 
transforms. Formally, these follow from those of the 7£-function. However, the Toda limit is not 
easy to control analytically, and the direct derivation of the eigenvalue equations is not too hard 
and quite illuminating. 

Our results in Section [5] have some overlap with earlier results by Kharchev, Lebedev and 
Semenov-Tian-Shansky [19], who obtained functions closely related to F T (r];x,y) from the 
viewpoint of harmonic analysis for Faddeev's modular double of a quantum group [5]. The 
nonrelativistic nonperiodic Toda eigenfunctions are widely known as Whittaker functions, and 
meanwhile it has become customary to call eigenfunctions for q-Toda Hamiltonians Whittaker 
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functions as well. In particular, q-Whittaker functions were introduced by Olshanetsky and 
Rogov for rank 1 (their work can be traced from [2D]) and by Etingof for arbitrary rank |21j . 
and these functions have been further studied in various later papers (see e.g. |22] and referen- 
ces given there). We would like to stress that these functions are quite different from the ones 
at issue here and in [19]. The crux is that the former are only well defined for q not on the 
unit circle, whereas here and in [19] the eigenfunctions have a symmetric dependence on two 
generically distinct q 7 s, given by 

q + = exp(i7ra + /a_), g_ = exp(i-7ra_/a + ). (1-19) 

This state of affairs is closely related to the different character of the trigonometric gamma 
function (more widely known as the g-gamma function, with the restriction \q\ ^ 1 being indis- 
pensable) and the hyperbolic gamma function (which depends on parameters a + and a_ in the 
right half plane) . 

In Section [6] we study the nonrelativistic limit of the representations of the relativistic eigen- 
function, arriving at two distinct representations for the nonperiodic Toda eigenfunction that 
have been known for a long time. Just as for the relativistic case, its property of being also an 
eigenfunction for a dual Hamiltonian seems not to have been observed before. (These duality 
features are the quantum counterparts of duality features of the pertinent action-angle maps, 
first pointed out in [17].) To control one of the two pertinent limits, a novel limit transition for 
the hyperbolic gamma function is needed, whose proof is relegated to Appendix [Bl 



2 The 7£-function as a special case of the H-function 

2.1 The functions 1Z and TZ r 

The i?-function (jl.ip is defined as a contour integral over a variable z, with the z-dependence of 
the integrand encoded in a product of eight hyperbolic gamma function factors. (See Appendix lAl 
for a review of the relevant features of the hyperbolic gamma function.) Specifically, with 
suitable restrictions on the eight variables, the R- function is given by 



R(c;v,v) = J J^F(c ;v, z)K(c; z)F(c ;v, z)dz. (2.1) 



Here we have 

c = (co + ci + c 2 + c 3 )/2, (2.2) 
t-i/ t G(z±y + id-ia) 

F (d;y,z ) = -^7T — r—j — r^-, (2.3) 
G[±y + id — ia) 

(with f(w ± y) denoting f(w + y)f(w — y)), and K is given by 
with new parameters 

si = c + ci - a_/2, s 2 = c + c 2 - a+/2, s 3 = co + c 3 . (2.5) 

Also, recall a and a are defined by (jl.4|) . 

We do not need the definition of the contour C for general variable choices (this is discussed 
in I and Section 4 of the survey [10]); instead we presently define C for the cases at issue. For 
the special c-choice in (jl.2p we can use the duplication formula (jA.lOp to obtain 

#((6,0,0,0);*) =£(&;*), (2.6) 
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where 

r( b . z \- 1 G(2ib-ia) G{z + ib-ia) 
{ ' ' ~ G{z + to) G(ib - ia) G(2z + 2ib-ia)' { ' ' 

Using also the reflection equation (|A.6[) we deduce that 7Z is given by 

„ /, /ax 1 / 2 G(2ib — ia) „, , , „, 

7£(6; x, y) = — -4— , - t& + ia)G{±y 2 - ib 2 + ia) 

\z7r/ G(^o — ia) 

G(z ±x + ib- ia)G(z + ib - ia)G(z ±y/2 + ib/2 - ia) ^ 

G{z + ia)G(2z + 2ib - ia) Z ' ^ ' ' 

For the variable choice that is most relevant for Hilbert space purposes, namely, 

a+,a-,b,x,y > 0, (2.9) 

the contour C may be chosen equal to the real line in the z-plane, indented downwards near so 
as to avoid a pole of IC(b; z). From ()A.17p -( jA.18p it follows that the poles of /C(6; z) are located 
on the imaginary axis at 

z-z kl = 0, z - z k i + ib = ia + /2, ia_/2, 0, fe,ZeN. (2.10) 
Thus they are above the contour, whereas the remaining z-poles of the integrand at 

z + z H = ±x-ib, ±y/2-ib/2, k,leN, (2.11) 
are below C. 

From the above representation it is immediate that 7Z is symmetric under the interchange of 
the parameters a+ and a_: 

1Z(a + ,a-,b;x,y) = 7£(a_, a+, b;x, y). (2-12) 

It is not at all clear, however, that 1Z is also symmetric under the interchange of the positions x 
and y: 

1Z(a + , a_, b; x, y) = TZ(a + , a_ , b; y, x). (2-13) 

This self-duality feature follows in particular from a second relation between 1Z and R, namely, 

K(b; x, y) = R{{b, b, 6, b)/2; x/2, y). (2.14) 

(This is equation (4.8) in [9].) Indeed, this second c-choice yields the same function fC(b; z) as 
the first one, so that substitution of (12. ip (with the same contour C) now yields (12.8P with x 
and y interchanged on the r.h.s. 

There are two more c-choices that lead from R to 1Z, namely, (6,0,6,0) and (6,6,0,0). 
Specifically, from equations (4.6) and (4.7) in [9] we have 

1Z(a + ,a-,b;x,y) = R(a+, 2a_ , (6, 0, 6, 0); x, y), (2.15) 
TZ(a+ , a_, 6; x, y) = R(2a-,a + , (6, 6, 0, 0); x, y). (2-16) 

From the definition of the i?-function we then obtain alternative integral representations for 
the 1Z- function from which the self-duality property (|2. 13|) is manifest. (Indeed, since we have 
c o = Co = 6 for these two choices, the integrand is invariant under the interchange of x and y.) 
On the other hand, the modular invariance property (I2.12p is not at all clear, since the integral 
representations involve the hyperbolic gamma function with a_ replaced by 2a_. Using (jA.lip . 
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they can be re-expressed in terms of the modular invariant function G(a + , a~;z). However, the 
resulting integrand is then still not modular invariant. Since it seems not to simplify and does 
not look illuminating, we do not detail it any further. 

The analyticity properties of the i?-function are known in great detail from Theorem 2.2 in I, 
cf. also Section 4 in the survey [ID]. Combining this theorem with the definition f)1.2|) of 1Z and 
its self-duality property (|2.13p . we deduce in particular that 1Z extends from the intervals (12.9f) 
to a meromorphic function in b, x and y, whose poles in x and y can only occur at the locations 

dzz = Ha — ib + Zki, z = x,y, k,l £ N. (2-17) 

We proceed to list further consequences of the i?-function theory for 7Z. Two features that 
are clear from each of the above integral representations are evenness and scale invariance (given 
scale invariance of G): 

K(b;x,y) =K{b;5x,5'y), 6,5' = +,-, (2.18) 
1Z(a + , a_, b; x, y) = lZ(\a + , Aa_ , Xb; Ax, Ay), A > 0. (2-19) 

A less obvious feature is the explicit evaluation 

K(b;x,ib) = 1. (2.20) 
It follows from the formula 

R{c;v,ic ) = 1, (2.21) 

(cf. equation (3.26) in I or Section 6 in |10j). by using any of the four relations (II. 2p . (12,14p . 
([2T5]) . ([2T6]) . Defining next 

TZ n (x) =TZ(b;x,y n ), y n = ib + ina-, n E N, (2.22) 

the eigenvalue AAE (analytic difference equation) for A + (b;y) entails 

s+{yn-ib) , v s + (y n + ib) 

— - — K n -i(x) H — - — K n+1 (x) = 2c + (x)TZ n (x). (2.23) 

s+iVn) s + (y n ) 

In view of (|2.20p . it follows from this that lZ n is of the form 

K n {x) = Pn{c+{x)), (2.24) 

where P n (z) is a polynomial in z of degree n and parity (— ) n . The relation of these polynomials 
to the g-Gegenbauer polynomials and to the Askey-Wilson polynomials associated with the four 
relevant c-choices is detailed at the end of Section 4 of [9j. 

The renormalized ^-function lZ r given by (jl.6p is the counterpart of the renormalized -R- 
function R r obtained from (|2.ip by omitting the product J} - G(isj) in K, cf. (|2.4p . (To see this, 
use (|Al3l) and (lATTOlU Clearly, it shares the features (^T2l) . ([fl5]l . (pT8l) and (pT9l) of the 
7?.-function, whereas (|2.20p is replaced by 

n r (b;x, i b)= ^~ ia \ . (2.25) 
G(2^o — ta) 

The renormalizing factor in the function lZ r ensures that it has no poles that are independent 
of x and y, cf. Theorem 2.2 in I. More precisely, 1Z r (a + , a_, b; x, y) extends to a function that is 
meromorphic in the domain 

D+ = {(a + ,a_,b,x,y) G C 5 | Rea+ > 0, Rea_ > 0}, (2.26) 

and whose poles can only occur at the locations (|2.17p . 
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It is not obvious, but true that for the special 6-choices b mn (|1.15p we have an equality 

7Zr(b mn ; x, y) = M(b mn ; x, y) + M(b mn ; -x, y), (2.27) 

where M(b mn ;x, y) is the function defined at the end of Section III in our paper [UJ. Therefore, 
lZ r (b;x,y) is the continuous (indeed, real-analytic) interpolation to arbitrary b £ K of the 
function given by equation (3.74) in [TT], which is defined only for the b- values b mn . (Note the 
latter are dense in M when the ratio a + /a„ is irrational.) 
It will not cause surprise that in the free case we have 

M (b 00 ; x, y) = exp(iaxy/2). (2.28) 

It would take us too far afield, however, to detail all of the functions M(b mn ;x,y), m,n G Z. 
For our purposes it is enough to specify their general structure: They are 'elementary' in the 
sense that they can be written 

M(b mn ;x,y) = exp(iaxy/2)R mn (e + (x), e_(x), e + (y), e_(y)), (2.29) 

where R mn is a rational function of its four arguments, cf. Section III in [TT]. In Subsection 12.41 
we deduce this structure in another way (namely, by exploiting parameter shifts). Moreover, 
for the case where m and n are not both positive or both non-positive, this structure can be 
understood from the novel Fourier transform representations (|3.48p - (|3.5ip . cf. Subsection 14.11 

We postpone the proof of the equality assertion (|2.27|) to Subsection l2.31 An ingredient of this 
proof is the asymptotic behavior of TZ r (b;x,y) as x goes to oo, and this is most easily obtained 
as a corollary of the asymptotics of a closely related function E(b;x,y), defined by (12.381) . 



2.2 The functions E and F 



The function E(b;x,y) can be viewed as a specialization of the function denoted £( , y;v,v) in II 
and [10]. The relation between 7 and c reads 

(2.30) 



7(a+, a_, c) = (c - a, c\ - a-/2, c 2 - a+/2, c 3 ). 
In particular, the 'free' case c = yields 
7f = (-a, -0-/2, -a+/2, 0). 



(2.31) 



The switch from c to 7 is crucial for uncovering further symmetries: The function £(pj\v,v) 
is invariant under transformations on 7 (i.e., permutations and even sign changes), cf. II. 
(In [B] this -D4 symmetry has been reobtained in a quite different way.) It is defined by 



£(r,v,v) 



-Rr(r,v,v). 



Here, the generalized (BC\) Harish-Chandra c- function is given by 

3 

- J 



(2.32) 



c(7; v) 



1 



G(2v + ia) 



(2.33) 



/Lt=0 



the dual of 7 by 

7 = Jj, 



( 1 1 
1 1 

1 -1 

Vi -1 



1 1 \ 

-1 -1 

1 -1 

•1 1 J 



(2.34) 
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and the constant by 

X(t) = exp (iapy • 7/4 - (a+ + at + a+a_)/8]), a = 27r/a+a_. (2.35) 
Denoting the 7-vectors corresponding to the two one- parameter families 

c = (6, 0,0,0), (6,6,6, 6)/2, (2.36) 
by 7 j "f^K it is easy to check that 

J 7 (1) =7 {2) - (2-37) 
Defining 

E(b;x,y)=£( 7 W;x,y/2), (2.38) 
a straightforward calculation (using the duplication formula (|A.10p ) yields 

E(6; x, y) = - ^ K r (b; x, y), (2.39) 
c{b;x)c{b;y) 

where we have introduced the constant 

<p{b) = exp{iab(b - 2a)/4), (2.40) 
and generalized (Ai) Harish-Chandra c- function 

G(z + ia — ib) . , , 

= G(z + u,) • (2 ' 41) 

Recalling (I2.14D and using (I2.37D . it readily follows that we also have 

E(b;x,y) =£( 7 M;x/2,y). (2.42) 

It involves more work to obtain the relations between E and £ corresponding to (|2.15p 
and (ETTBJ) . Setting 

7^) = 7 (a+,2a_, (6,0,6,0)), = 7 (2o_ > o +> (6,6,0,0)), (2.43) 
these are given by 

E(a + ,a_,b;x,y) = £(a + , 2a_, 7 (3) ; x, y) = £(2a_, a+, 7 (4) ; x, y) . (2.44) 

(These formulas amount to special cases of the doubling identity for the ^-function obtained in 
Section 6 of [9].) 

The relation (|2.39p between E and 7Z r yields a similarity transformation turning the AAOs 
(|1.13p into A±(b;x), A±(b;y), where 

As(b; z) = c(6; z)~ l A s (b; z)c(b; z) = I£_ 4 + V s (b; z)Tl ia _ s , 6 = +,-, (2.45) 

V 5 (b; z) ee S6(z + ib)ss(z-ib + ia- S )^ 
ss{z)ss{z + ia-s) 

From this it is easy to verify that these AAOs are formally self-adjoint operators on L 2 (M) (by 
contrast to the AAOs (|1.13p ). and that they are invariant under the transformation 

6 (->■ a+ + a_ - 6. (2.47) 
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It is not obvious, but true that we also have 

E(b;x,y) = E(a+ + a_ -b;x,y). (2.48) 

This symmetry property can be derived from (|2.38[) and the -D4 invariance of the ^-function: 
We have 

7 ( 1 ) = (6-a,-a_/2,-a+/2,0), (2.49) 

so the map b 1— > 2a — b yields a 7-vector related to 7W by a sign flip of the first and last 
component. 

Combining (|2,39p and (|2.4ip with the analyticity features of 1Z (cf. the paragraph con- 
taining (|2.17]) ) ; we deduce that E(b;x,y) is meromorphic in b, x and y, with ^-independent 
pole locations 

z = — 2ia — Zki, z = x,y, fc, I £ N, (2.50) 

corresponding to the factor G(x + ia)G(y + ia), and 6-dependent poles at 

z = ib + Zki, z = —ib + 2ia + Zki, z = x,y, k,l£N. (2.51) 

The main disadvantage of the function E(6; x, y) compared to the 7£ r -function is that it is 
not even in x and y, since the c-functions in (|2.39p are not even. Instead, it satisfies 

E(6; -x, y) = -u(b; x)E(b; x, y), (2.52) 

where 

u(b; z) = -c(b; z)/c(b; -z) = -G(z ± (ia - ib))/G{z ± ia). (2.53) 

On the other hand, E inherits all other important properties of lZ r , and is the simplest function 
to use for Hilbert space purposes. In particular, it has the 'unitary asymptotics' 

E(6; x, y) ~ exp(iaxy/2) — u(b; —y) exp(—iaxy/2), (2.54) 
y£ (0,oo), x — > 00, 

cf. Theorem 1.2 in II. Here, the ^-function encodes the scattering associated with the AAOs 
A±(b;x), reinterpreted as commuting self-adjoint operators on the Hilbert space 
More precisely, using corresponding results on the ^-function from II and III, it follows that the 
generalized Fourier transform 

F : C = C^({0,oo))cH = L 2 ({0,oo),dy)^H = L 2 ((0,oo),dx), (2.55) 

defined by 

/ a \ i/ 2 f°° 

(jT^Xx) = (— J E(b;x,y^(y)dy, i> G C, (2.56) 

extends to a unitary operator, provided the coupling b is suitably restricted. Specifically, it 
suffices to require 

&€[0,a+ + a_]. (2.57) 

The self-adjointness of the operators A±(b) on H associated to the AAOs A±(b; x) for b in this 
interval can then be easily understood from the unitarity of J-: they are the pullbacks to "H 
under T of the self-adjoint operators of multiplication by 2c±(y) on 'H. 
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As already mentioned, these statements follow from II and III by specialization, but it may 
help to look first at Section 9 in the survey [TDJ. Starting from the representation (|2.38p . the 
vector belongs to the polytope P given by equation (9.2) in [ID], provided 6 G (0,2a). 
Therefore, the transform associated with £(j^; v, v) (defined by equation (9.4)) is an isometry. 
As a consequence, the transform (|2.56p is an isometry. (The normalization factor in (|2.56p differs 
from that in equation (9.4) in |10] to accommodate the scale factor 1/2 in the y-dependence 
of 8 in (|2.38p .) Isometry of the inverse transform is then clear from the self-duality of E(6; x, y). 
Next, for 6 = we have the identity 

E(0; x, y) = K r (0; x,y) = 2 cos(axy/2), (2.58) 

cf. (|239jl - (rX4Tj l and equation (7.33) in [10] (also, note that for b = we have 7W = 7f, cf. (TXIDD 
and (|2.3ip ) . In view of the symmetry (|2.48|) , it follows that J- amounts to the cosine transform 
for 6 = and b = 2a, so these transforms are unitary as well. More generally, we obtain a family 
of unitary operators 

T(a + ,a-.,b), (a+, a_, b) £ U u = (0, oo) 2 x [0, a+ + a_], (2.59) 

which is strongly continuous on the parameter set H u and satisfies 

J 7 (a + , a_, 6) = J r (a + , a_ , a + + a_ — 6), (2.60) 

cf. Theorem 3.3 in III. 

The G-function asymptotics ()A.13P entails that the c-function (|2.41j) has asymptotics 

c(6; x) ~ <f){b) ±l exp(=Fa6x/2), Re (x) -)■ ±00, (2.61) 
with 4>(b) given by (|2.40p . Hence the u-function (|2.53j) has asymptotics 

u(b;x) <P(b) ±2 , Re (a) -> ±00. (2.62) 

Also, the reflection equation (1A.6[) and the complex conjugation relation (]A.9[) entail 

u(6;-s)«(6;x) = 1, |u(6;x)| = l, 6,xGM. (2.63) 
Thus, if we set 

F(6; x, y) = 0(6)- 1 (-u(6; x))V2(_ u (6 ; y)) 1 /2 E (6 ; x , y ), 6, x , y > 0, (2.64) 

(with the square root phase factors reducing to 1 for 6 = 0), then F has asymptotics 

F(6;x,y) ~ [— u(b; y)] 1 ^ 2 exp(iaxy/2) + [— u(6; y}~ 1 ^ 2 exp(— iaxy/2), x — > 00, (2.65) 

and if we replace E by F in the above unitary transform (I2.56p . we retain unitarity. 
Introducing the weight function 

w(b;x) = l/c(6;±x) = G(±x + ia)/G(±x + ia - ib), (2.66) 

we have 

w(b;x)>0, 6,x>0, (2.67) 
and we can also write F in terms of lZ r as 

F(6; x, y) = w(b; x) l / 2 w(b; y) l/2 TZ r (b; x, y), 6, x, y > 0, (2.68) 
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with the positive square roots understood. Hence, F is a joint eigenfunction of the AAOs 
H±(b;x) and H±(b;y) with eigenvalues (|1.14|) . where 

H s (b; z) = w(b; zf^Agib; z)w(b; z)' 1 ' 2 

From the G-AAEs (1A.2|) we deduce 

w(b; x) = 4s + (x)s_(x)w r (b; x), w r (b; x) = G(±x — ia + ib). (2.70) 

For b G (0,2a) the reduced weight function w r (b;x) is positive for all real x, and since it is 
also even, its positive square root for x > has a real-analytic extension to an even positive 
function on all of K. By contrast, it is clear from (|2.70p that wfax) 1 / 2 , x > 0, extends to an 
odd real-analytic function on M. As a consequence, one can also view the transform associated 
with F(b; x, y), b € (0, 2a), as a unitary transform from the odd subspace of L 2 (M, dy) onto the 
odd subspace of L 2 (]R, dx). This is the viewpoint taken in [12], where we studied this transform 
(among other ones) for the special 6- values Na+ with JV £ N*. As shown there, for b > 2a 
unitarity and self-adjointness generically break down in a way that can be understood in great 
detail. 

To be sure, the precise connection between the above functions F and lZ r and the functions F r 
and E r from [12] is not clear at face value. But the latter are derived from the functions 
M((N + l)a + ;x,y) of [IT] , as specified below equation (1.42) in [12], so this connection is 
encoded in the identities (|2.27p for the special cases (m,n) = (N + 1, 0), N € N. 



2.3 The identities (12.271) and their consequences 

We proceed to prove the general identities (|2.27j) . Our reasoning involves in particular a com- 
parison of the behavior for x — > oo of the functions on the l.h.s. and r.h.s. For 7Z r (b; x,y) this 
asymptotics easily follows upon combining fj2.39[) . (|2.54p and (|2.6ip : 

lZ r (b; x, y) ~ exp(— abx/2) c(b; ry) ex.p(riaxy /2), (2-71) 

T=+- 

b e K, y > 0, x -> oo. 

Next, we consider the functions M(b mn ; ±x, y). To begin with, they are eigenfunctions 
of the four AAOs ([1713]) (where b 

— bmn) with eigenvalues (j!.14j) . cf. Theorem II. 3 in [11] . 
Their 'elementary' form (|2.29p follows from equations (3.65)-(3.68) in [11]. The function 
Kn+,N- («+ 5 O- ; y) occurring in these formulas is specified in equation (3.2), with Sn s given by 
equation (2.21). In turn, the coefficients in equation (2.21) are defined via equations (2.2)-(2.5) 
in [11]. (See also Subsection 14.11 for more information on these special cases.) It is straightfor- 
ward to obtain the asymptotics for Re (x) — > oo from these explicit formulas. Specifically, this 
yields 

M(b mn ; ±x, y) = e W (-»bmnx/2)c(b mn ; ±y)e ±iax y' 2 [l + 0(e^ Rc <*>)] , (2.72) 
y > 0, Re (x) — > oo. 

The decay rate p is the minimum of the two numbers 27r/a±, and the implied constant can be 
chosen uniform for Im (x) varying over M. 

Comparing (I2.7ip and (I2.72p . it follows that the functions on the l.h.s. and r.h.s. of (12.271) 
have the same asymptotics for x — > oo. It therefore suffices to prove that for fixed a + , a_, y > 
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and m, n G Z, they must be proportional as functions of x. Moreover, we may as well assume 
a+/a_ is irrational, since equality for this case entails equality for all a+,a_ > 0. (Indeed, 
the functions M(b mn ; ±x,y) are manifestly real-analytic in a + and a_ for a+,a_ > 0, and this 
real-analyticity property is also valid for lZ r (b; x,y), cf. I.) 

The key consequence of the irrationality assumption is that the vector space of meromorphic 
joint solutions f(x) to the AAEs 

A±(ma + + na_; x)f(x) = 2c±(y)f(x), m^nGX, y > 0, a + /a_ ^ Q, (2.73) 

is two-dimensional. To explain why this is so, we first note that the functions M(b mn ; ±x, y) 
are independent solutions to (12.73p . their independence already being clear from their general 
form ()2.29p . Moreover, it follows from their uniform asymptotics ()2.72p that there exists a po- 
sitive number A, depending on the fixed variables a+, a_, m, n and y, but not on x, such that 
in the half plane Re (x) > A both functions are zero-free, and satisfy 

lim M(b mn ;x,y)/M(b mn ;-x,y) = 0, Re (x) G (c_, c+) C [A, oo). (2.74) 

Im (x)— >oo 

We are now in the position to invoke a result from Section 1 in [23], to the effect that the 
above suffices for any joint meromorphic solution f(x) of (|2.73p to be a linear combination of 
the two functions M (b mn ; ±x , y) . Since TZ r (b mn ; x,y) is an even meromorphic joint solution, 
the functions on the l.h.s. and r.h.s. of (|2.27p are proportional, so their equality now follows. 
In particular, for the free case m = n = we recover the identity (12.58P from (|2.27p - (|2.28p . 
Moreover, taking y = ib mn in (|2.27p . we can invoke f|2.25|) to deduce the corollary 

•, x , r /, •, \ G(i(m-l/2)a + +i(n-l/2)a-) , n nr .. 

M(b mn ;x,ib mn ) + M(b mn ;-x,ib mn ) = — ^ / + 77^^- 2 " 75 

G(i(zm — l/2)a + + i{2n — l/2ja_) 

Using the G-AAEs (fO|) . the r.h.s. can be rewritten in terms of sine- functions. For the special 
case m = N + 1, n = 0, the resulting identity amounts to equation (2.78) in |11| . cf. also 
Subsection 14.11 

We would like to add in passing that it is very plausible that (|2.74p is not necessary for two- 
dimensionality. Indeed, denoting by V c the field of meromorphic functions with period c G C*, 
any third independent joint meromorphic solution would have to be both of the form 

f(x) = pi(x)M(b mn ;x,y) + p 2 (x)M(b mn ; -x,y), Pi,P2&V ia+ , (2.76) 

and of the form 

f(x) = qi(x)M(b mn ;x,y) + q 2 (x)M(b mn ; -x,y), q\, q 2 G V{ a _. (2.77) 

Since the intersection of the fields V% a+ and Vi a _ reduces to the constants when a+/a_ is 
irrational, we expect (but are unable to prove) that this simultaneous representation should 
lead to a contradiction without appealing to (|2.74p . 

Now that we have proved ()2.27p . it follows that the function 1Z r (a+, a_, b; x, y), which is 
real-analytic on the parameter set 

II = {(a+, d-,b) G (0, oo) 2 x M}, (2.78) 

is the continuous interpolation of the functions on the r.h.s. of (|2.27p . which are only defined 
for the dense subset of 'elementary' parameters 



II e i = {(a + , a_, b) G II | b = ma + + na-, m, n G Z}. 



(2.79) 
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The natural question whether another linear combination of M(b mn ; ±x, y) that is independent 
from the even one admits a continuous interpolation as well remains open. In this connection we 
should point out that our reasoning at the end of Section 3 of [23] renders this extremely unlikely, 
but is not conclusive. Indeed, we cannot rule out that the sequence of functions Q- given by 
equation (3.15) in [24] . with N + € N, gives rise to an infinity of distinct limits L_, corresponding 
to distinct subsequences. (This oversight is of no consequence for the later sections in [24] .) For 
the same reason, the analogous assertion about the R- function, made at the end of [10] . has not 
been completely proved. 

Before turning to parameter shifts, we derive a non-obvious reality feature of lZ r from the 
relations (|2.27p . namely, 

7Z r (a + ,a_,b;x,y)=7Z r (a + ,a-,b;x,y), V (a + , a_, b, x, y) £ II x 1R 2 . (2.80) 

The point is that from the explicit formulas for the functions M it is apparent that we have 

M(b mn ;x,y) = M(b mn ;-x,y), x,y £R, (2.81) 

cf. equation (3.73) in [11]. Therefore, (|2.80p is clear from (|2.27p and interpolation. As a corollary, 
this yields reality of 1Z and F for real parameters and variables, cf. (|1.6p and (I2.68p . (Alter- 
natively, this reality property of the 7£-function follows from that of the iZ-function proved in 
Lemma 2.1 of III.) 

2.4 Parameter shifts 

We continue to summarize results concerning parameter shifts from [8], inasmuch as they apply 
to the present A\ context. In Section 1 of [8] we introduced the up-shifts 

^ u) (-) = ^y(^-^_J, (2-82) 

satisfying 

sf{x)A &l {b-x) =A 5 ,{b + a s -x)sf\x), (2.83) 
and the down-shifts 

Sf (b; x) = -A- [s s (x - ib)s s (x + ia. s - ib)Tf a ,-(*-> -i)] , (2.84) 
s s (x) 

satisfying 

Sf ) (b;x)A s ,(b;x) = A s ,(b - a s ;x)S^ d) (b;x), (2.85) 

where 5,5' = +,— . Clearly, the up-shifts S^"\x) and S_ (x) commute, and the down-shifts 
S+\bi;x) and S^fo'^x) commute as well. The shifts are also related by 

Si u) (x)sf(b;x) = A 5 (b;x) 2 - 4co S 2 (7r(b - a-s)/a s ), (2.86) 
S { s d \b + a s ;x)S ( s u \x) = A s (b;x) 2 - 4 cos 2 (irb/a s ), (2.87) 

where 5 = +, — . It is a matter of straightforward calculations to verify the formulas f|2.83|) 
and (gjjgD - (gI57D . 

Starting from the joint eigenfunctions exp(±iaxy/2) of ^4±(0; x) with eigenvalues 2c±(y), one 
can now obtain joint eigenfunctions with the same eigenvalues for A±(b mn ; x) by acting with 
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the shifts on the plane waves. By construction, these joint eigenfunctions are of the elemen- 
tary form (|2.29|) . Choosing a + /a_ irrational, it follows from two-dimensionality of the joint 
eigenspace that these eigenfunctions are (generally y-dependent) multiples of M(b mn ; ±x;y). 
A more telling action of the shifts is encoded in 

Sf\x)1l r {b-x,y) =4ss{y + ib)s s (y-ib)K r {b + a s ;x,y), (2.88) 
sf\b-x)n r {b;x,y)=K r {b-a 5 ;x,y). (2.89) 

Indeed, these relations hold for arbitrary b. For b = b mn , it then follows by using ()2.27p that 
they also hold for the summands M(b mn ; ±x, y). (This is because their plane wave factors are 
independent, cf. QZjggJ) .) 

The equations (|2.89p and f|2.88|) follow from a suitable specialization of equations (3.11) 
and (3.13) in [5]. But in the present A\ case we can also derive them quite easily by using 
the elementary joint eigenfunctions M(b mn ;x, y) with a + /a_ ^ Q. Indeed, once we have shown 
that (|2.89p . (|2.88p hold for y > 0, b = b mn , m,n G Z, and with 1Z r replaced by M, it is easy 
to deduce (|2.89p . ([2.88P from (|2.27p and interpolation. (Note that the four shifts commute with 
parity.) Their validity for these special cases can be readily verified: One need only show that 
the functions on the l.h.s. and r.h.s. have the same x — > oo asymptotics, and using (|2.72p this 
causes little difficulty. 

Next, we obtain the counterparts of the Ag(b; x)- and 7£ r -shifts for the AAOs As(b; x ) and 
their joint eigenfunction E(b;x,y). (For the BC\ setting we did this in Section 8 of [9]; as in 
previous cases, it is in fact simpler and more illuminating to obtain the relevant formulas by 
direct means, instead of by specialization.) They are given by 

s ^ l)s 7d^) sf)(l)c[bil) - (2M) 

Sf(b:x) ^ — ^ -Sf> (6; *)«;(&;*). (2.91) 

c(b — as; x) 

A moment's thought shows that this entails the validity of (|2.85p - ([2.87p with S, A replaced 
by S, A. Also, using the definition (l2lTj) of the c-function and the G-AAEs (fAT2l) . we obtain 
the explicit formulas 



4"<k*) -is,- f-.'ff;'*^ (2 - 92) 

s s (x)s s {x + ids) 



c W (h . rr\ - t x s s (x + ib)s s (x + ib - ia-g) 

S s (b,x)-T ia _ s ss{x)ss{x + ia _ s) -ia-s ■ ( 2 -93) 

Notice that they imply 

si u) (2a-b-x)=S { s d \b-x). (2.94) 
Finally, a straightforward calculation yields the counterparts of (|2,88p and (|2,89D : 

sl u) {b;x)E(b;x,y) = 2e s (-iTTb)s 5 (y + ib)E{b + a 5 ;x,y), (2.95) 
S^ d \b;x)E(b;x,y) = 2es(iir{b - a_ s ))ss{y ~ ib + ia_s)E(b - a s ;x,y). (2.96) 

To conclude this section, we point out that the eight shifts acting on x have duals acting 
on y given by the formulas QgJjlD , d^Mjl and <^Uh - H2~m\) with x -)• y. By self-duality, their 
respective actions on lZ r (b;x,y) and E(b;x,y) follow from the above by interchanging x and y. 
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3 Five minimal representations of the ^-function 

We begin this section by focusing on the kernel function 

K(b;x,v) :=G({±x±v-ib)/2) = G(($iz + 5 2 v - ib)/2). (3.1) 

Si,S 2 =+- 

We have established that this function satisfies three independent kernel identities. We expect 
that these might be useful in other contexts than the present one. Indeed, here we only need the 
special case (|3.14p of the first of the identities. We collect the three identities in the following 
proposition. 

Proposition 3.1. Letting b, d £ C, we have 
ss(x — ib + id) 



s 5 {x) 

s$(v — id) 
ss(v) 

ss(x — ib) 



K(b; x — ia_s, «) + (*—>■ —i) 

K(b;x,v — ia_s) + (i—>-—i), (3-2) 



s s (x) 

_ ss(v - ib) 



ss(v) 
8 S ((x - ib)/2) 



K(b; x — 2ia_$, «) + («'—>■ —i) 

K{b; x,v- 2ia- 5 ) + (t -»• -i), (3.3) 



ss(x/2) 

_ s 5 ({v-ib)/2) 



K(b; x — ia_$, v) + (i — > —i) 

K{b\x,v-ia-8) + {i-*-i). (3.4) 



ss{v/2) 

Proof. To prove (|3.2p . we divide l.h.s. and r.h.s. by K(b;x — zo_^,^) and use the AAEs (|A.2|) 
to write the result as 

sg(x — ib + id) ^ ss{x + ib — id) cg((x + v — ib)/2) cg((x — v — ib)/2) 
sg(x) ss(x) cs((x + v + ib)/2) cg((x - v + ib)/2) 

_ s 5 (v - id) c s ((x - v - ib)/2) | s s {v + id) c s ({x + v- ib)/2) 
s s (v) c 5 ((x - v + ib)/2) s s (v) c s ((x + v + ib)/2) ' 

Both sides are 2m,5-periodic functions of x with equal limits 

ej(±(— ib + id)) + es(±(— ib — id)), Rex — > dboo. (3-6) 

The residues at the (generically simple) poles x = 0, x = ia$ in the period strip clearly cancel. 
By Liouville's theorem, it remains to check that the residues at the poles x = ±v — ib/2 ± ia$ 
cancel as well, and this is a routine calculation. 

Next, we divide (|3.3[) by K (b; x, v) and use (1A.2|) to obtain 

sg(x — ib) c$((x — ias/2 ± v + ib)/2) 

+ (x — > —x) 



s$(x) c$((x — ia^gj2 ± v — ib)/2) 

ss(v — ib) cs((v — ia_&j2 ± x + ib)/2) 



+ («-)•-«). (3.7) 



ss(v) cs((v — ia_s/2 ± x — ib)/2) 
Both sides are 2m,5-periodic functions of x with equal limits 

rj--u\ . (s- 'u\ s s{v~ib) s s (v + ib) 

e 5 (±ib) + e ( 5(=F^0) = 1 — , Rex-^ioo. (3.8) 

ss{v) s 5 (v) 
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The residues at x = and x = ias manifestly cancel. It is a straightforward calculation to verify 
that the residues at the remaining poles x = ±ia$ ± ias/ 2 ± v + ib cancel, too. Hence (|3.3|) 
follows. 

Finally, to prove (|3.4p . we divide both sides by K(b;x — ia_s,v) and use (|A.2|) to get as the 
counterpart of (|3.5[) : 

s 5 ((x-ib)/2) ssjjx + tb)/2) c s {{x + v- ib)/2) c 5 {{x -v- ib)/2) 
s s (x/2) s s {x/2) c 5 ((x + v + ib)/2)c 5 ((x-v + ib)/2) 

= s s ((v-ib)/2)c s ((x-v-ib)/2) s s ((v + ib)/2)c s ((x + v-ib)/2) 

s s (v/2) c s ((x-v + ib)/2) s s {v/2) c 5 {{x + v + ib)/2) ' [ ' ' 

Both sides are 2m^-periodic functions of x with equal limits 

e s (±(-ib/2)) + e s (±(-3ib/2)), Rex^ioo. (3.10) 

As before, residue cancellation at x = and x = ias is immediate, whereas the verification 
that the residues at the remaining poles x = ±ia$ ± v — ib/2 cancel as well involves a bit more 
work. ■ 

From (ll.7p we see that the identity (|3.2p can be rewritten as 

As(b- d;x)K(b;x,v) = A s (d;v)K(b;x,v), 5 = +,-. (3.11) 

At first sight, one might think that the identities (|3.3p and (|3.4p can also be rewritten by using 
a rescaled version of the two commuting A\ difference operators A±(b;z). The two difference 
operators 

s s ((z-ib)/2) s 6 ((z + ib)/2) 

ss(z/2) ^-' + s 5 (z/2) ' — ' ' (3 - 12) 

featuring in (|3.4p . do not even commute, however. Thus no such rescaling is possible for (|3.4p . 
The difference operators 

ss{z~ib) ssjz + ib) 

corresponding to (|3.3p do commute. Even so, one can only rescale one of the operators such 
that it takes the A\ form (jl.7p . but not both at once. 

For the purpose of studying the A\ operators, then, we can only make use of (|3.2p . More 
specifically, our starting point is the special case d = 0: 

A 5 (b;x)K(b;x,v) = + H ia jK{b; x,v). (3.14) 

In order to exploit this identity, we introduce the Fourier transform 

B(b;x,y) = - I dvK(b;x,v)exp(iavy/2), b G (0,2a), x,y €M.. (3.15) 

The integral is well defined, since the 6-restriction ensures that the f-poles of the integrand at 

±v = x + 2ia — ib + zm, ±v = —x + 2ia — ib + z^i, k, I £ N, (3.16) 

(cf. (|A,17p -( |A,16P ). stay away from the real axis, and since the G-asymptotics (|A.13P entails an 
exponential decay 

K(b;x,v) ~ exp(=Fa6i)/2), b > 0, x G C, Ret;-^±oo. (3.17) 
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These features also imply that B extends from the real x-axis to a function that is holomorphic 
in the strip Imx G (—2a + b, 2a — b). 
Next, we temporarily assume 

6g(0, a s /2), a s = min(a+, o_). (3.18) 

Then the action of the shifts in the AAOs A±{b\ x) given by (|1.7p is well defined on B(b; x, y), 
provided we restrict x to a strip |Imx| < a s /2. Moreover, we may take the shifts under the 
integral sign in ()3. 15j) and use the kernel identity ()3.14|) to obtain 

A$(b; x)B(b; x, y) = - j dv 'S~\ K(b; x, v + ria_g) exp(i7rm//a + a_), (3.19) 

|Imx| < a s /2. 

Upon shifting contours K — > M. ± ia_£, no poles are met, and so we obtain 

- es(ry) / dvK(b; x, v ) exp(iirvy/a+a_). (3.20) 

2 T =+- JR-ira^s 

The integrands of both terms are now equal, and the contours can be shifted back to M without 
changing the value of the integrals. Hence we deduce the eigenvalue equations 

A±(b;x)B(b;x,y) = 2c±{y)B(b; x,y), \lmx\ < a s /2. (3.21) 

Reverting to our previous assumption b G (0,2a), we proceed to obtain two different repre- 
sentations of B(b; x, y). To this end we use the Plancherel relation 

/ dpf(p)g(p) = f dqf(q)g(-q), f,g€ L 2 (R) D L\R), (3.22) 

with the Fourier transform defined by 

h(q) = / dpexp(iapq)h(p), h = f,g. (3.23) 
Jm. 

Rewriting B as 

B(b,x,y)-^dp G{p _ [x _ lh)/2)G{p+{x + ib)/2) > ^-24) 
b G (0,2a), x,yeR, 
we now define 

G(p+(x-ib)/2) G(p + (x-ib)/2) 

h{p) = G[p-{x-ibW Hp) = G(p + (x + ib)/2y (3 ' 25) 

/ \ _ G(p - (x + ib)/2) exp(iapy) _ G{p - (x + ib)/2) exp{iapy) 
91{P) = G(p + (x + l b)/2) ' 92{P) = G(p-(x-ib)/2) • (3 - 26) 

We can calculate the Fourier transforms of these four functions by using Proposition lC.il Doing 
so, we use the Plancherel relation (|3.22p and then replace q by z + y/2 to obtain the two 
representations announced above: 

x ^ ,x f , G(z±(x-y)/2-ia + ib/2) . ^ 

B{b;x,y) = G{±x + ia-ib) / dz ) "t, ^. JJr, 3.27 

J R G(z ±(x + y)/2 + ia- ib/2) 
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B(b; x, y) = G{ia - ib) 2 I dzG{±z ± y/2 - ia + ib/2) exp{iazx). (3.28) 

Jr 

Next, we compare (|3.28p to (|3.24p . deducing that it can be rewritten as 

B(b; x, y) = G(ia - ib) 2 B{2a - b; y, x). (3.29) 
Also, defining a new function C(b;x,y) by (jl.lip . we see that (|3.27p amounts to 

/27T 

BCb; x, y) = J —G(±x + ia - ib)C(2a - b; x, y). (3.30) 
V a 



The function C(a+, a_, b; x, y) is of central importance for what follows. Writing it as 

faT I" yr G(a + ,a^;z - u s ) 
2ti" Jr A, G(a + ,a_;z - d s ) 

where we have introduced 

u ± = ±(x-y)/2 + ib/2, d± = ±(x + y)/2 - ib/2, (3.32) 

we infer that its behavior under analytic continuation in its 5 variables is immediate from the 
general analysis in Appendix B of I (with iV specialized to 2). We proceed to summarize the 
salient information. 

To this end, we need the function i?(a+,a_;z) discussed in Appendix El cf. (IA.20p . Intro- 
ducing 

P(a + , a_, b; x, y) = C(a+, a_, b; x, y)E(a+, a_; ±x + ib — ia)E(a + ,a_; ±y + ib — ia), (3. 33) 

the product function P(a + ,a-,b;x,y) extends from (0,oo) 2 x (0, a+ + a_) x R 2 to a function 
that is holomorphic in the domain 

D(a+,a-,b) = {(a + ,a-,b,x,y) G C 5 | Rea + > 0, Rea_ > 0, Re(6/a + a_) > 0}. (3.34) 

Hence C is meromorphic in Z)(a+,a_,6), with poles occurring solely at the zeros 

±x = 2ia — ib + zm, ±y = 2ia — ib + z^i, k, I G N, (3.35) 

of the -E-product, cf. ()A.2ip : moreover, the maximal multiplicity of a pole at z = zq, with 
z = x,y, is given by the zero multiplicity at z = zq of the pertinent E-factor. 

The corresponding meromorphy properties of B are now clear from its relation (|3.30p to C: 
It continues meromorphically to the domain D{a+, a_, 2a — b). From (|3.29p we then deduce that 
B{aj r ,a-,b;x,y) has a meromorphic extension to the larger domain D + (I2.26p . Using (13.301) 
again, it now follows that C has a meromorphic extension to D + as well. 

We proceed to obtain further information on the function C. First, we list features that are 
immediate from its definition (II. lip and properties of the G- function, cf. Appendix (Aj 

C{a + ,a_,b;x,y) = C(a-,a + ,b;x,y), (modular invariance), (3.36) 
C{a + ,a-,b;x,y) = C{a + ,a-,b;y,x), (self-duality), (3.37) 
C(b;x,y) = C(b;5x,5'y), 5,5' = +,—, (evenness), (3.38) 
C(a + , a_, 6; x, y) = C{\a + , Aa_, Xb; Xx, Xy), X > 0, (scale invariance), (3.39) 
C(a + , a~, b; x, y) G R, o + ,a_>0, b G (0,2a), x,yGR, (real- valuedness). (3.40) 
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Clearly, the relations (|3.36|) f)3.39[) are well defined and hold true on D + (|2.26p . Also, the 
property (|3.4U|) can be rendered manifest by substituting the integral representation ()A.5|) in 
the four G-functions and combining factors to obtain 

, , / a f , ft dw smixw) sm(yw) cosh(6w) s\n(2zw) \ .„ 

C(a + ,a.,b;x,y) = J — / dz cos / V . ; f \ ■ u \ " ( 3 - 41 

V 27T y K yjjj w smn(a+u;) smn(a„wj / 



x exp 



dw f cos(xw) cos(yw) sinh(ou>) cos(2zw) b 



l jo \ sinh(a + u;) sinh(a_u?) a + a_w 



where (a + ,a_,b,x,y) G (0, oo) 2 x (0,2a) x M 2 . 

Second, combining (|3.37p with (13. 29ft and (13.30p . we deduce 

G(ia- i b)C(b;x,y) _ G(ib - ia)C(2a - b;x,y) (6 . symmetry) . (3 . 42) 



G(x + ia — ib)G{y + ia — ib) G(x — ia + ib)G(y — ia + ib) 

Third, we can use Proposition IC . 1 1 once more to obtain from (jl.lip the explicit result 

C(b;x,ib) = G(ia — 2ib)G(ib — ia) 2 , (normalization), (3.43) 

where (a + ,o_,6, x) G (0, oo) 2 x (0,a) x R. Last but not least, we claim that we have the joint 
eigenvalue equations 

A ± (b;x)C(b;x,y) = 2c ± (y)C(b; x,y), A ± (b;y)C(b; x,y) = 2c ± (x)C(b; x,y). (3.44) 

To prove this claim, we first note that the eigenvalue equations (|3.21|) continue meromorphi- 
cally to D + . Next, we observe that the G-AAEs (|A.2p imply 

G(±x + ia - iby 1 A 8 {b;x)G(±x + ia - ib) = A 5 (2a - b;x), 5 = +,-. (3.45) 

(Here, we view the l.h.s. as the product of three operators acting on meromorphic functions.) 
Hence we obtain via (|3.30p 

A±(2a-b;x)C{2a-b;x,y) = 2c±{y)C{2a - b; x,y), (3.46) 

which is equivalent to the first two AAEs in (|3.44[) . The last two are then clear from the 
self-duality relation (13.371) . 

All of the properties of C just derived also hold true for the function G(ib — ia)lZ r defined 
by (II. 6p and (ll.2p , cf . Section [2l By using solely the eigenvalue properties (|3.44p , the evenness 
properties (|3.38[) . and the normalization (|3.43p . we can now show that these two functions 
coincide, as announced in the Introduction, cf. (|1.12p . Specifically, applying the uniqueness 
argument in Subsection 12.31 to C in its dependence on x, we obtain the equality ()1.12p up 
to a proportionality factor p(a + ,a-,b,y). Repeating this argument for the y-dependence, we 
see that the proportionality factor can only depend on the parameters a + , a_, b. From the 
normalization relation (|3.43p it then follows that p = 1, thus proving (11 . 12|) . 

Let us now collect the resulting minimal representations of the 7£-function. From ()1.12p 
and (11.61) we obtain 



nib- x y) = 1 G (™-™) [ d G(z±( X -y)/2-ib/2) 

X > V) ^a-Z G{ib - ia? k G(z ±(x + y)/2 + ib/2) ' ^ A() 

Next, combining (|3.24[) and (|3.30p . we deduce 

. 1 G(2ib — ia) . 

TZ(b; x, y) = —— — jG{±x + ia - ib) 

y/a + a- G(ib — ia) A 
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f G(z±x/2- ia + l b/2) 

J R ° G(z ± x/2 + ia- ib/2) ex W az y>' V A *> 
and using (|3.29[) we infer 

K(b\x,y)= -j G(2ib - ia)G(±x + ia - ib) f dz °}f ^ V ' 2 ~ 2 expfiazx). (3.49) 

Vo+o^ Jr G(z ± y/2 + 

Finally, using the self-duality property of 7£, we obtain from (|3.48p and (|3.49p the representations 



1 G(2ib - ia) 
yja + a- G(ib — ia)' 



K(b; x, y) = __ 7 A- 7 -— -±G(±y + ia - ib) 



G(z±y/2-ia + l b/2) 
d ~ G(z ± y/2 + ia - ib/2) e W a ^' ^ bU ) 

K(b; x, y) = 1 G(2ib - ia)G{±y + ia - ib) 

V a + a - 

, G(z± x/2 -ib/2) , , 
dz G{z±x/2 + l b/2) ^azy). (3.51) 

The five representations (|3.47p - (|3.5ip are well defined and hold true for (a + ,a^,b,x,y) G 
(0,oo) 2 x (0,2a) x M 2 . (Combining ([3^7]) with (13^2]) . we get a further representation that 
we do not consider.) 

Taking stock of the above developments, we note that we might have started from the first 
minimal representation (|3.47p to define the ^-function. Then many of its properties follow 
quite easily. On the other hand, it seems not feasible to give a direct proof of its crucial joint 
eigenfunction property. With hindsight, however, this can be shown by first obtaining the 
second representation (|3.48p (say) via Proposition IC.ll and then using the identity (|3.14p to 
arrive at (I3.2ip . From this the joint eigenfunction property (I3.44h follows as before. 

Another important property of 1Z is its asymptotic behavior for x — > oo. Like other features 
addressed in this section, this can already be gleaned from Section [21 via the specialization of 
the more general asymptotics of the function £(7; v, v) obtained in Theorem 1.2 of II. However, 
provided we restrict b to the interval (0, 2a), it is quite easy to obtain the x — > 00 asymptotics 
directly from the new representations of 7Z in terms of a Fourier transform. 

To detail this, let us first note that we need only consider the function E(6;x,y), which 
we can now view as being defined via (|2.39p - (l2.4ip . (Indeed, there is no difficulty in obtai- 
ning the asymptotics of the c- function; in this connection, compare (|2.6ip . (|2.54p and (|2.7ip .) 
Using (I3.49p . we deduce that for 6 6 (0, 2a) the E- function has the representation 

E(b- x y) - *@ G{lb ~ m) ° {x + m) [ dz G(z±y/2-ib/2) 
l ' ' V) ~ ^aTa-Z c(b;y) G(x - ia + ib) J M Z G(z ± y/2 + ib/2) ex ^ iazx >- ^ bl > 

Letting y € (0,oo), we can shift the contour up by a — 6/2 + e, where e > is small enough so 
that only the simple poles at 

z = ±y/2 - ib/2 + ia, (3.53) 
are encountered. The residues at these poles easily follow from (|A.19p . yielding a contribution 

M(6; x)(exp(iaxy /2) + c(6; — y) exp(— iaxy /2)/c(b; y)), (3.54) 
with the multiplier given by 



M(b; x) = <t>{b)G{x + ia)G(-x + ia - ib) exp(-ax(a - 6/2)). 



(3.55) 
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Now from (|A.13P we see that M(b;x) converges to 1 for x — > oo. To recover the asymp- 
totics (|2,54|) . therefore, it is enough to show that the r.h.s. of (|3.52|) with z replaced by 
z + ia — ib/2 + ie vanishes for x — > oo. 

To prove this, we write the shifted contour integral as 

1 G(ib-ia) f G(z + ie ± y/2 - ib + ia) 



f G(z + ie ± y /2 — ib + ia) 

M(b\ x) expl — eax) / dz — — ; — exp (lazx . 3.56 

Jr G(z + ie±y/2 + ia) Fy 1 y ' 



Now one need only use (IA.13j) to verify that the integrand is bounded by a multiple of 
exp(— ab\z\/2), which implies the function (|3.56[) does converge to for x — > oo. 

We stress that this short argument only yields (|2.54p under the restriction b € (0,2a). In 
particular, by contrast to the previous contour integral representation used in II, one must cope 
with an inevitable contour pinching when one tries to use (|3.52p to go beyond this 6-interval. 

Another issue is that stronger asymptotic estimates than just obtained are necessary to 
recover the Hilbert space transform features for the E-function sketched in Subsection 12.31 
cf. (I2.55p - (l2.60p . It is beyond our scope to study this further, but we would like to repeat that 
the 6-interval [0, 2a] cannot be enlarged without losing the critical unitarity and self-adjointness 
properties [12]. 

At face value, the new representations (|3.47p - (|3.5ip seem to hold promise for a direct proof 
of the shift properties of 1Z r , cf. (I2.88p - (|2.89p . Even so, we were unable to push this through. 
To date, therefore, the only reasoning yielding the properties for general b is to first derive them 
for special b- values, as sketched in Subsection I2.41 We now turn to a study of the 7£ r -function 
for these special values. 



4 Specializations and nonrelativistic limit 

4.1 Elementary special cases 

As already mentioned, the functions 

R N (a + ,a-;x,y) = Kr(a + ,a-,(N + l)a + ;x,y), N e N, (4.1) 

have been extensively studied before. They were first obtained more than twenty years ago [25], 
and then reconsidered from an algebraic and function-theoretic viewpoint in [11] and from 
a representation-theoretic viewpoint in a paper by van Diejen and Kirillov [26]. The corre- 
sponding Hilbert space transforms were studied in great detail in \12\ . 

Our first goal in this section is to demonstrate how the elementary character of these functions 
can be directly understood from the Fourier transform representations (|3.48p - (|3.5ip . Indeed, 
thus far the relation encoded in (I2.27P has only been shown by appealing to a uniqueness 
argument. The crux is that for the choices b = (N + l)a± one can use the G-AAEs (|A.2p to 
obtain integrals that can be explicitly evaluated by a residue calculation. 

Specifically, let us start from ()3.49p to obtain first 

4 -7V-l 

R N (x, y) = G(i(N + l)a+ - ia)G(±x + ia - i{N + l)a+) 

V a + a - 
r N 

x / dze iazx /T[c-(z±y/2-i(N - 2j)a+/2). (4.2) 
Jr j=0 



Now we recall that (|3.49|) is valid for b G (0,2a), which implies we have Na+ < a_ in (|4.2p . 
Taking y > from now on, it follows that the integrand has 2N + 2 simple poles in the strip 
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Imz 6 (0, a_). The product yields a function that is ia_-periodic in z. Thus, denoting the 
integral by In, we have 

N ( 1 \ 
/at - exp(-27rx/a + )/Ar = 2vri ^ ^ Res ( z = - {by + ia_ + i(N - 2j)a+) j . (4.3) 

5=+ - j=o ^ ' 
The residues are easily calculated, and hence we obtain 



(_t)"+ia. 
s+(x) 



( e i«yx/2 N e_(-(jV-2j» + _ ^ 

, s -(v) ~L EI s^(y+i(k-j)a + )sin(ir(j - k)a+/a^) 



(4.4) 



Also, the prefactor can be calculated by using (|A.2j) once more, combined with ()A.12p . Intro- 
ducing 

N 

P N (z)= _ 2s-(z + ija + ), z = x,y, (4.5) 
j=-JV 

we get 

9-2AT-1 i \ N 

^frn 2sin ^' a +/ a -)- ( 4 - 6 ) 

a_ P N (x) AX 

For iV = the product of (fOj) and yields 
o / \ sin(7rxj//a + a-) 

More generally, the product can be written as 

R N (x,y) = (-i) N+1 (K N (x,y) - K N (x,-y))/P N (x)P N (y), (4.8) 
where we have set 

v 



Kx(x,y) = exp(ixy/a + a_) J^2sin(7rZa + /a_) 



1=1 

N N 
N e -(( 2 J ~ N ) x ) EI s-(y-ika + ) ]J S-(y + ika + ) 

xV fc= ^. , n k f-> +1 . (4.9) 

7^ EI sin(7r(j - A;)a+/a_) 

Introducing the phase factor 

g = exp(i7ra + /a_), (4-10) 

it is not hard to see from (j4.9|) that if at is of the form 

K N (x,y) = eyip(ixy/a + a^)e-(Nx + Ny)S N (q;e-(-2x),e-(-2y)), (4.11) 
N 

S N (q;r,t)= J2c^\q)rH l , (4.12) 



fc,«=0 



with cij (g) a rational function of (/. 
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Thus far, our conclusions about Rn and Kn were only based on an explicit evaluation of 
the representation f|3.49|) . However, a lot more information follows upon using the features of 
lZ r (b;x,y). In particular, the function K^ix^y) / Pn{x)Pn{v) is a joint eigenfunction of the four 
AAOs (|1.13p (where b = (N + l)a+) with eigenvalues (|1.14p . since this holds true for Rn(x, y). 
Also, the self-duality and evenness of i?Ar imply 

K N (x, y) = K N (y, x) = K N (-x, -y), (4.13) 

and this entails that the coefficients in (|4,12p have the symmetry properties 

C (N) _ (N) _ (N) k , _ N (4 U) 

c kl ~ Hk ~ c N-k,N-V K,l — U, . . . ,1\ . K^- 1 ^) 

Of course, this can be easily checked for small N, but for arbitrary integers (|4,14p is not at all 
obvious from (|4.9p . 

One more feature of the coefficients is that they are Laurent polynomials in q with integer 
coefficients. Like the symmetries (|4.14|) . it is not a routine matter to show this from (|4.9|) . 
The crux is, however, that the above functions Kn coincide with those of [IT], by virtue of 
the uniqueness argument used in Subsection 12. 31 The coefficients were studied in detail in 
Section II of [TTJ, and there the interested reader can find explicit formulas for the coefficients 
as Laurent polynomials in q. See also the paper by van Diejen and Kirillov [26], where yet 
different representations of the functions K]y(x,y) were obtained. 

With a little more effort, the elementary character of lZ r for the more general b- values 

6+_ = (N + l)a+-Ma-, JV.MeN, (4.15) 

can also be understood from (|3.49p . Indeed, from (IA.2P it follows by a straightforward calculation 
that we have an identity 

M 

Yl2c + (v + i((N + l)a + + (M + l-2k)a-)) 
G(v - ib+-/2) = fe= i 

G{v + ib + _/2) n , „ [ ' 

V ' ' EI 2c_ {v + Mo. + {N — 2j)a+)) 

3=0 

Using this identity several times (together with (|A.12p for the factor G(ib^ — ia)), we deduce 

from (|3.49p and (jl.6p the representation 

N M 

Y[ 2sin(j7ra + /a_) 2s + (x — zfca_) 

K r (b + -;x,y) = -L^ *^JT ( 4 - 17 ) 

] 2 sin(A;7ra_/a + ) 2s_(x — ija + ) 

k=i j=~N 

M 

f] 4c+ (z+l(±y + i(N + l)a + + i(M + l-2k)a-)) 

dze 1 ' 



jiazx k—1 



N 

TT 4c_ (z + i( ± y + iMa_ + i(N - 2j)a + )) 

3=0 

The denominator of the integrand has no zero for z£i unless M is odd and N is even. The zeros 
of the corresponding factor s-(z ±y/2) are then matched by the zeros of the factor s+(z ± y/2) 
of the numerator. After a suitable contour shift, we can expand the numerator product into 

exponentials, yielding a sum of convergent integrals (recall we require b^ £ (0,2a)). When M 

is even or ./V is odd, we can do the same without a contour shift. Each of the integrals is then 
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basically of the same form as previously evaluated for the case M = 0. (More in detail, the same 
2N + 2 poles arise in the period strip for the c_-product.) 

From these observations the general structure anticipated in (|2.27j) and (|2.29p readily follows, 
provided n > and m < 0, or m > and n < 0. In Section III of [11] we studied the 
functions ()4.17|) in considerable detail, but it is beyond our scope to derive the explicit form 
used there directly from their representation (I4.17p . We do add that it seems plausible that 
the factorization exhibited in equations (3.3)-(3.4) of |llj can be understood by a more refined 
analysis of the above sum of contour integrals. In any case, we repeat that equality of the 
pertinent functions follows from the uniqueness argument explained in Subsection 12.31 



4.2 The nonrelativistic limit 

We begin this subsection with a remark addressed to physicist readers, who may care about 
dimension issues. In our paper |llj . which we had occasion to cite several times in the previous 
subsection, the variable y of the present paper was denoted by p. This is a widely used notation 
for the spectral variable in nonrelativistic quantum mechanics, where p is viewed as a momentum. 
In our relativistic setting, however, it is far more natural to view the scale parameters o+ and a_ 
as having dimension [position], and then the 'geometric' and 'spectral' variables x and y both 
have dimension [position] as well. (To be more specific, one of the scale parameters can be 
viewed as an interaction range, and the other one as the Compton wave length h/mc of the 
relativistic particles under consideration.) This goes to explain our change from p to y. 

Of course, from a mathematical viewpoint such notational issues and the notion of dimension 
may be ignored. When taking the nonrelativistic limit, however, these physical considerations 
point the way. We need to let the speed of light c go to oo, so one of the scale parameters should 
go to 0. In particular, we cannot retain modular invariance. Accordingly, we first set 

a+ = 2-ir/fi, a_ = h/3, {j,, h,fi>0. (4.18) 

Here, we view /3 as 1/mc, with m the particle mass, and we trade a+ for a parameter fi with 
dimension [position] -1 to avoid a great many factors tt. Next, the spectral variable y (dual 
position) is replaced by the momentum variable 

P = M/P. (4.19) 

Finally, the coupling parameter b (with dimension [position]) is replaced by 

b = g(3, g>0, (4.20) 

so that g has dimension [action] = [position] x [momentum] as well as Planck's constant h. 
With these changes in effect, it is easy to verify the expansion 

A+(b;x) = 2 + p 2 A + 0(^), /3->0, (4.21) 

where 

A := -h 2 d 2 x - ghficoth(fix/2)d x - c/V/4, (4.22) 
and the limits 

lim A„ (b; x) = exp(-mg/h)T^ + (i -> -i) =: M, x > 0, (4.23) 

limA + (6; y ) = + (i -> -*) =: i, (4.24) 
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while A_(b;y) has no sensible limit. The eigenvalue of A + (b;x) on lZ(b;x,y) satisfies 

2c + (y) = 2 + /3V/4 + 0(/3 4 ), (4.25) 

so that the eigenvalue of A becomes p 2 /4, while the eigenvalues of the monodromy operator M 
and dual AAO A remain 2 cosh.(irp / hfi) and 2cosh(^x/2), resp. 

Likewise, the similarity-transformed AAOs A± (|2.45p and Hamiltonians H± (|2.69p yield 



A+(b;x),H + (b;x) = 2 + /3 2 F + 0(/3 4 ) , -»■ 0, (4.26) 
4sinh 2 (/xc/2) 



« : = + ,T,2/^ - (1-27) 



lim^-(6;x) = limff-(6 ;;c ) = T5, /( ,+Tf 2j , /(l =:>l, I > 0, (4.28) 
lim^; ,) = T* w + Eim^Ezm ... A, (4.29) 

= E_ (™f^ (^f == A- (4.30) 

(We take this opportunity to point out that the L BC\ version' of the monodromy operator (14.281) 
specified in equation (8.20) of jlO] contains an error: the phase factor on the r.h.s. should be 
replaced by 1.) 

The operators A, M and H, Ai are related by a similarity transformation with the limit 
function 

lim <f>(b)/c(b;x) = w nr {g/h; fix/2) 1 / 2 , x > 0, (4.31) 

/3->0 

where 

w nv {X;r) = (2sinhr) 2 \ A,r>0, (4.32) 
the operators A and A by similarity with 

lim G(ia-2ib)G(ib-ia)/c(b;y) = 2/c nT (g/h;p/hfi), (4.33) 

/3->0 

where the Harish-Chandra c-function is given by 

and the operators A and H by similarity with w nr (g / h; p/ hfi) 1 ^ 2 , where 

w ni (X;k) = l/c nr (X;±k). (4.35) 

The limits (|4.3ip and (|4.33p can be readily verified via the definition (|2.4ip of the relativistic 
c-function and the above reparametrizations (14.18p - (|4.20p by using the G-limits ()A.25[) ( 1A.26|1 . 
The functions (|4.32p and (|4.34p are normalized so that 

lim w nr (\; r) = 1, lim c nr (A; k) = 1. (4.36) 
A— >o A— >o 

Next, we study the nonrelativistic limit of the 7£-function. To this end we first use the 
definition (jl.ip . the scaling property (|2.19|) and the limit 

lim R(tt, t, tc; v, tu) = 2.F1 (cq + m, Co — iu, cq + C2 + 1/2; — sh 2 v) , (4.37) 
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c = (co + ci + c 2 + c 3 )/2, (4.38) 
established and discussed in I. We write the limit 

limTl(2ir/ijL,H/3,gP;x,Pp/fJ,) =: ip nr (g/h; fix / 2, p/ hp), (4.39) 

/9-)-0 

in terms of the dimensionless quantities 

X = g/h, r = fix/2, k = p/hfi, (4.40) 
already used in (14.321) and (|4.34j) . The result reads 

^ nr (A; r, k) = 2 F 1 ((X + ik)/2, (A - ik)/2, A + 1/2; - sinh 2 r) . (4.41) 

(See [27] for a limit that is related to cf. (j!T24) .l 

Likewise, the alternative representation (|2.14|) entails 

V> nr (A; r, fc) = 2 iq (A + ik,X- ik, A + 1/2; - sinh 2 (r/2)) , (4.42) 

whereas (|2.15p - (|2.16p again give rise to (|4.41 [) . The equality of (|4.4ip and (|4.42|) can be rewritten 

as 

2 F 1 (a,b,a + b+ l/2;4iu(l - w)) = 2 Fi(2a,2b,a + b + 1/2; w), (4.43) 

which is a well-known quadratic transformation, cf. e.g. [28^ p. 125]. 

Using other familiar features of the hyper geometric function, it is not difficult to verify that 
the operators A, M and A do have the expected eigenvalues p 2 /4, 2 cosh(p/hfi) and 2 cosh(//a;/2) 
on the limit function ip m (g/h; jjLx/2,p/hp). More specifically, for A this amounts to the ODE 
satisfied by 2 F\ and for A this involves the contiguous relations. The M-eigenfunction property 
follows by using the known analytic continuation of 2 F\(a, b, c; w) across the logarithmic branch 
cut w G [1, oo). 

The above specialization of the hyper geometric function basically yields the so-called conical 
(or Mehler) function. To be specific, the latter can be defined by 

^^(coshr) = 2A ^r^Vl/2j 2Fl(A + iK X ~ iK X + V2; (1 " coshr )/ 2 )' ( 444 ) 

cf. |16t equation (14.3.15)] and the hypergeometric function occurring here equals the one 
in fljggp . 

We now consider the nonrelativistic limit of the minimal representations of TZ(b; x, y) derived 
in Section [3J We were not able to obtain a sensible limit for the second one, given by (|3.48D , 
For the remaining four, however, the limit can be handled in a sense to be explained shortly. 
For expository reasons we first list the resulting representations for ip nr (X; r, k): 

t it _J__„( k ^(^±mX\, (4.45) 



2 A r(A) 2 J R (coshr + cosht) A \ \cosh({t - r)/2) 
r(2A)(sinhr) 1 ~ 2A f T({it - \±ik + c + l)/2 
2 2X + 1 ir J R t T({it + X±ik + c+l)/2 

r(2A) 

4vrr(A) 2 r(A ± ik) 



dt-^- — ' , " ; exp((i(t-ic)r), c> X — l, (4.46) 



/ dtV{{it + A±tfc)/2)r((-tt + X±ik)/2)exp(itr), (4.47) 

f du ^tk) (448) 



2 x T(X±ik) J R (cosh r + cosh t) : 
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We proceed to discuss these formulas. First, we note that they are derived under the as- 
sumption 

A, r, k e (0,oo), (4.49) 

and that this implies that the integrals in (|4.45p . (|4.47p and f|4.48|) are absolutely convergent. 
The integral in (|4.46|) . however, is only absolutely convergent for A > 1/2; For A £ (0,1/2] it 
should be viewed as a Fourier transform in the sense of tempered distributions. 

Second, we compare these formulas to results in [16], where a host of representations for 2-F1 
and its conical function specialization are listed. Formula (|4.48p can be readily found there: It 
can be obtained from equation (14.12.4), which can be written 

P l/2 ~ X (cosh r) - , / ^Xsmhr)*- 1 / 2 /■» cos kt 

^-1/2^°^)-^ T{X±ik) J dt (cosh r + cosh t) x- ^-50) 

(This involves the duplication formula of the gamma function, equation (5.5) in [16] ■) As 
they stand, the three remaining representations do not occur in [16]. However, as was pointed 
out by a referee, they can also be tied in with results in the vast literature connected to the 
hypergeometric function. To begin with, formula (|4.47p can be derived (with some effort) by 
combining equations (15.8.14) and (15.6.7) in [16]. It seems that the formulas (14.45P and (14.461) 
cannot be obtained by using only [16] or some other standard reference book. Even so, they agree 
with known results. Indeed, (|4.45p amounts to equation (2.3) in the paper [29], whereas (14.461) 
can be derived by combining several sources. Specifically, the integral in (|4.46j) can be viewed 
as a special case of the Meyer G- function, cf. Section 16.17 in [16] and p. 144, (2) of [30]. After 
contour deformation, a residue calculation leads to a formula involving a linear combination of 
two 's with gamma function coefficients, cf. equations (16.17.2) and (16.17.3) in [16] or (7) 
in [30] . Finally, it follows from equation (3.2(27)) in |31] that the latter formula yields the 
conical function as represented by (|4.50p . 

Third, none of the four representations (I4.45p ~( f4.48p has been obtained with complete rigor. 
The difficulty is to obtain uniform tail bounds on the pertinent integrands that allow an ap- 
plication of the dominated convergence theorem. (In fact, to date a similar snag has not yet 
been obviated for the limit transition (|4.37p either.) In this connection we should add that we 
were unable to verify directly that each of the four formulas gives rise to a joint eigenfunction 
of the operators A, M and A with eigenvalues p 2 /4, 2 cosh(p / hfj.) and 2 cosh(/ia;/2). On the 
other hand, a direct check of the joint eigenfunction properties of the five relativistic represen- 
tations (|3.47p - (|3.5ip seems not feasible either. 

We continue by sketching the derivation of the four formulas (|4.45p - (|4.48p . First, we observe 
that any factor of the form 

G(a + ,a-]ia — ita~), (4-51) 

with t not depending on a_, can be treated via (|A.25p . Indeed, a scaling by a + yields 

G(l, k; i/2 + K (i/2 - it)) ~ — 2 ^ k y 0. (4.52) 

T{t) exp(i m(2iTK)) 

In particular, this yields not only the asymptotics of the numerical prefactors, viz., 
1 r(2A)exp(2Aln(/3fyi)) 



yfa+aZG{ia — 2ib) 27rf3h 
G(ia - ib) 2 /j,T(2X) 



P 0, (4.53) 



y/a^aZGiia - 2ib) T(\y 



P -± 0, (4.54) 
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but also that of the y-dependent prefactor in ()3.50[) (|3.51|) : 

nc -i. a. \ 27r/3%exp(-2Aln(/3n/i)) 

G^ct^&iy) r(A±zfc) , ft -)• 0. (4.55) 

To handle the x-dependent prefactor in (j3.48j) — (|3.49j) . however, (|A.25P is of no help and (|A.26|) 
seems not to apply either. But in fact we can use the G-AAEs (|AT2j) to first write 

G(x + ia- ib) _ 2 x Gj> - ia + /2 + ia_/2 - ib) ^ 
G(x — ia + ib) G(x — ia+ /2 — ia_ /2 + ife) ' 

and then (|A.26P can be invoked to deduce the dominant asymptotics 



G(±x + ia-i6)~exp — (2sinhr) i -^ A , x > 0, ft ->■ 0. (4.57) 

The plane waves in (|3.48p - (|3.51|) become 

exp(iazx) = exp(i/izx/ fth), exp(iazy) = exp(izp/H). (4.58) 

For the first and second plane wave we now switch to a new variable t given by 

z -> /3^/2, z -> (4.59) 

so that they turn into exp(itr) and exp(itk), resp. Likewise, in (|3.47p we change z to t/f/, to get 
dimensionless G-function arguments. 

With these variable changes in place, we proceed to look at the asymptotic behavior of the 
G-ratios in (I3.47p - (|3,5ip . For the first case this is immediate from (lA.26p . and this easily 
yields (|4.45p when the pointwise limit is interchanged with the integration. (As alluded to 
above, an L 1 -bound uniform for /3 near is needed to make the interchange rigorous. As 
well as in the next cases, no such bound is available for now.) For the second case (|3.48p . it 
seems not possible to get from the pointwise behavior of the integrand (with x > 0) as f3 goes 
to a factor exp(— nx//3ti) that takes care of the diverging factor exp(irx/ f3h) coming from the 
prefactor (I4.57p . By contrast, for the third case we may and will make a shift 

t-H + i/(2nPh)-ic, (4.60) 

with cGR chosen so as to stay away from poles while shifting and taking ft to 0. This cancels 
the diverging factor and results in (|4.46p via (|A.25P . Finally, an application of (1A.25P and (1A.26|) 
leads to the limits (|4.47|) and (j4.48j) . resp. 

To conclude this subsection, we point out that in view of (12.391) and (ll.2p the nonrelativistic 
limit of E(6; x, y) is given by 

2t^(A;r) 1 /2 ^ = . 

c nr (A; k) 

cf. (|4.3ip ~ (l4.34p . It has the unitary asymptotics 

E nr (A; r, k) ~ exp(ir/c) — u nr (X; —k) exp(-irk), r — > oo, (4.62) 
with the scattering function 

. n ,v c nr (X;k) T(ik)T(X-ik) . . 

Uni{X ' k) = "MA^y = ~T(-ik)T(\ + ik)' (463) 
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The latter is normalized so that it equals 1 for A = 1, just as u(b;z) f|2.53|) is normalized to 
equal 1 for b = a±. In this connection we would like to add that from (|4.7|) and (|2.39p - (|2.4ip 
one readily deduces 

E(a±;x,y) = 2i sin(7rxy/o + a_). (4.64) 

Hence the reparametrizations (j4.18fl — (|4.2Q[) and (|4.40p yield 

E m (l;r, k) = 2ismkr. (4.65) 

Accordingly, the 'free' theory with which the scattering is compared is given by the sine transform 
(and not by the cosine transform, which arises for b = A = 0, cf. f)2 .581) 1 . 

5 The relativistic Toda case 
5.1 Taking the relativistic Toda limit 

Throughout this section, we require that the parameters a+ and a_ be positive. It is also 
convenient to require 

(x,y) G (0,oo) 2 , (5.1) 

until further notice. In keeping with our outline in the Introduction, we begin by considering 
the b- values (|1.16p . In this case the w- function (|2.66p reads 

w(a - ij;z) = G(±z + ia)/G(±z - 7), a = (a+ + a_)/2, 7 e R, (5.2) 

and hence is no longer real-valued for real z and 7 / 0. By contrast, the ^-function (|2.53p is 
given by 

u(a-ir,z) = -G(z±'y)/G(z±ia), (5.3) 
so it is still unitary for real z; moreover, it is even in 7. The Hamiltonians (I2.69P can be written 

u ( ■ \ {cs(z + ia~s/2±j)\ 1/2 „ u , .> , . 

H s (a - i-y; z) = — ^ — - — — exp(ia_ s d z ) + (t ->■ -«), (5.4) 

V s$(z)ss{z + ia-s) ) 

so they remain formally self-adjoint for real z; they are also even in 7. 

Next, we consider the five representations of the joint eigenfunction F(a — i-y; x, y) of the four 
Hamiltonians H±(a — ij;x) and H±(a — i"f,y)- Combining (|2.68p and (|1.6p with (|3.47|) — (|3.5ip . 
these are given by 



G(- 7 ) ( G(±x + ia)G(±y + ia) y z f G(z± (x - y)/2 - ia/2 - 7/2) 

{ G(±x - 7) G(±y - 7) J h * G(z ± (x + y)/2 + ia/2 + 7 /2) ' 1 ' J 

G(- 7 ) /G(±a? + io) G(±y + m)\ 1/2 f , „, . ; . , 

vM= ( c U +7 ) C (4- 7) ) /„ ± 1/2 - • a ' 2 + ^ ^' a ^ p- 6 » 

( g<^ + ")^ + ") ) " 2 /• ^(i, ± y/a _ io/2 - 7 /2) e*p( ia zs), (5.7) 
^/alaZ V G(±x + 7) G(±y - 7) / ./r 

G(-7) (G(±x + ia) G(±y + ia)\ 1/2 f , , , , , 

/ ^ - \ , . A , . / dzG (d=z ± i//2 - io/2 + 7/2 exp iazx , 5.8 

G( 7 ) /G(±x + ia) G(±« + ia) \ 1/2 /" . , , . 

- UJ , ,. . , dzG(±z±x/2-ia/2->y/2)exp(iazy). 5.9 
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(The square roots are positive for 7 = 0.) As they stand, none of these representations yields 
a manifestly real- valued function for 7 / 0. However, comparing (|5.7|) and (|5.8|) . we see that 
these formulae are related by a complex conjugation (take z — > —z in one of them to check this). 
Likewise, (|5.6p and (|5.9p are related by a complex conjugation. Since the five formulae yield the 
same function F(a — 27; x, y), this function is in fact real-valued. 

This reality feature can be tied to the b — > 2a — b symmetry of the E- function, cf. (|2,48p . In- 
deed, the latter invariance implies that E(a — i'j; x, y) is even in 7. Now since the ^-function (|5.3p 
and the phase <p{a — vy) (given by (|2.40p ) are manifestly even in 7, it follows that F(a — ij; x, y) 
is even in 7, cf. (|2.64p . Comparing once again (|5.7p with (j5.8j) . and (|5.6p with ()5.9|) . we see that 
these formulae are also related by flipping the sign of 7, in accord with evenness. 

Substituting 

7-M? + A, x-tx + A, (5.10) 
we are now prepared to study the Toda limit A — > 00. First, for r£lwe have 
r c s (x + ira^ s /2 + rj + 2A)c^(x + ira_ s /2 - rj) . 

hm — ; = 1 + es{-2x - iras + 2r]). (5.11) 

A->oo ss(x + A)ss{x + ira_s + A) 

Thus we obtain relativistic nonperiodic Toda Hamiltonians given by 

Hj (77; x) = lim Hs(a — irj — iA; x + A) 

A— >oo 

= [1 + e s {-2x - ia^ 5 + 2r?)] 1/2 exp(ia- S d x ) + (i -i). (5.12) 

(The square roots are positive for x — > 00.) Clearly, these are formally self-adjoint on L 2 (R, dx). 
In this connection we point out that in view of the diverging x-shift, we may and will from now 
on allow x to vary over R in the Toda quantities, whereas we continue to require that y be 
positive. 

Next, we note 

lim e s (-2A)c s {y + ira_ 5 /2 + r, + A)c 5 (y + ira_ 5 /2 - rj - A) = e 5 (27?)/4. (5.13) 

A— s>oo 

Hence we get dual Hamiltonians 

H J(W V) = , lim e s (-A)H 5 {a - ii] - zA; y) 

= ^s S (yr 1/2 {eMia-sd y ) + ex V (-ia_ s d y )) ssiy)^ 2 , (5.14) 

which are formally positive on L 2 ((0, 00), dy). 

To obtain the Toda limit of the joint eigenfunction F(a — 27; x, y) involves a greater effort. 
The key tool is the asymptotics (|A.13P of the hyperbolic gamma function. This enables us to 
show that the limit 

F T (r]]x,y)= lim F(a — ir\ — iA; x + A, y) (5.15) 

A— >oo 

exists for each of the five representations f)5.5[) — f)5.9[) . The details now follow. 

To start with, the asymptotic behavior for A — > 00 of the five prefactors can be assembled 
from the three formulae 

G(x + v + A) ~ e~ ix exp + A) 2 + x 2 + 2x( V + A)) V (5.16) 

G(±y - V - A) ~ e 2 ^ exp (fa + A) 2 + y 2 )) , (5.17) 
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G(±(x + A) + ia) ~ exp(7raa(x + A)). (5.18) 

Next, consider the integrand of (|5,5p with the substitutions (|5.10p . Two of the four G- functions 
are invariant, and the remaining two yield 

G(z - (a - y)/2 - ia/2 - ij/2 - A) 
G(z + (x + y)/2 + ia/2 + rj/2 + A) 

^ e 2^ exp ^([ z + y /2]2 + [( x + ia + r? )/ 2 + A] 2 )^ . (5.19) 



If we now combine the A-dependent terms coming from the prefactor in (|5.5|) . then we see that 
they cancel the A-dependent term in (|5.19p . The product of the remaining terms is readily 
verified to be given by 

° i±y + m) ^ V2 exp(3, X /2) - U ! '[ 2 ^ ~ exp (%m) , (5.20) 



a + a.G(x-r])J * K 1 G(z - (x + y)/2 + ia/2 + rj/2) " \ 4 
with 

M = 2z 2 + 2zy - y 2 /2 - iax - a 2 /2 + iarj. (5.21) 
Finally, since we integrate z in ()5.20|) over R, we may shift z by y/2, yielding the limit function 

1/2 



F (jj;x,y) 



G{± ^ m \ ) exp(3z X /2) exp (y 2 - ia (x - r,) - a 2 /2)) 
J dzG{±z + (x-r)- ia) /2) exp (^-(z 2 + 2zy) ^ . (5.22) 



Turning to the integrand of (|5.6p . the substitution (|5.10p again leaves two of the four G- 
functions unchanged, as well as the plane wave factor. The remaining G-product has asymptotics 

G(±z + (x — ia + V )/2 + A)~e- 2i Xexp(-^(z 2 + [(x-ia + V )/2 + A] 2 )^ . (5.23) 

Combining this with the asymptotics of the prefactor following from (I5.16p - ()5.18p . the A- 
dependent terms drop out. Taking z — > — z in the resulting limit function yields 

1 /2 

F^; x, y) = ( gfe±jj ) exp(-3f X /2) exp ( ^ (y 2 + ia(x - r?) - a 2 /2)) 
\a + a^G(-x + rj) J \ 4 v '/ 

x y dzG(±z- (x - i] + ia) /2) exp (^^-{z 2 + 2zy)^j . (5.24) 

Comparing this representation to (|5.22p . we see that each of the factors on the right-hand side 
is matched by its complex-conjugate. Thus, the equality of f|5.22|) and (|5.24p is in keeping with 
the real-valuedness of F T (r];x,y) following from its being the limit of a real-valued function. 

Proceeding in the same way for (|5.9p . we obtain as its limit again (|5.22p . Of course this 
should be expected, since the factors on the right-hand side of (|5.9p and (|5.6p are related by 
complex conjugation. On the other hand, the equality of the limits of (|5.5p and (|5.9p yields 
a nontrivial check of the substantial limit calculations. 

At face value, the representations (|5.7p and (|5.8p seem not to give rise to a sensible Toda 
limit. In fact, however, they do, but it is expedient to postpone the details. First, we rewrite the 
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representations (|5.22p and (|5.24p in a more telling form, by bringing in the G-cousins Gl and Gr, 
cf. (|A.27p ~ (|A.28p . Indeed, a straightforward calculation yields the equivalent representations 

F T (??; x, y) = ( + \ V2 e ^ 2 /4 f dzGR (± z + {x _ ^2 _ ia /2)e iaz v, (5.25) 

and 

F r (7/;s,j/) = ( G ^ + ZQ ) V 2 e -^/4 /" cfoG^-^-^-ia^e^. (5.26) 
\a + a-G L {ri - x) J J R 

(Here, the square roots are positive for x — > oo, cf. ()A.3ip .) 

Now we use once again the Plancherel relation, as encoded in (|3.22j) ()3.23j) . Taking first 

f(p) = G R (p + (x- tj)/2 - ia/2) exp(ia[p + (x - rj + ia)/2) [y/2 + ia/4]), (5.27) 
g(p) = G R (-p +{x- rj)/2 - ia/2) exp(ia[-p + (x - r) + ia)/2] [-y/2 + ia/4]), (5.28) 

we can use ()C.45P with s = a/4 to compute the Fourier transforms. This yields 

1 /2 

/(?)= (?) e " W4_Mx ^( Z ^?( !B -^ + io )) G! i(9 + I'/ 2 - 3i o/4), (5-29) 
g( q ) = f—\ 1 e -W4-2ix exp f— q ( x - v + ia)) G L (-q - y/2 - 3m/4). (5.30) 



a J \2 
The integral in (|5.25p is therefore equal to 



exp f-^-(ia)(x - rj + iaU ^- f(q)g(-q)dq 

/ dqexp (-ia(q + ia/4) (a; - 7? + ia))G L (q ± y/2 - 3m/4). (5.31) 



e -i7r/2-4ix 



The new representation thus obtained can be somewhat simplified by reverting to the G- 
function, and by shifting the contour down by a/4 (recall y > 0). In this way we obtain 
from (I5.25P an alternative representation 

1 /2 

F T (r?; x, y) = f + ) eX p (za y 2 /8) 

\a + a~G r{x - ri) J 

x / dzG(z±y/2-ia)exp(-ia[z(x-r ? ) + z 2 / 2 ])- (5-32) 

./R+iO 



Proceeding in the same way for (|5.26p . we arrive at a fourth representation, namely 

x / dzG(z±y/2-ia)exp(-ia[z(x-r ? )-z 2 / 2 ])- (5-33) 

./R+iO 

Comparing it to (|5.32[) . we deduce once again real-valuedness of the function F T (rj;x,y) on 
R 2 x (0,oo). 

Having these two new representations at hand, we can see with hindsight that they can also 
be obtained from (|5.8|) and (|5.7|) . respectively. Indeed, when we let 

z + ia/2-r]/2- A/2 (5.34) 

in the integrand of (|5.8p . then we obtain the two G- functions featuring in (|5.32|) . times two 
A-dependent ones. If we now proceed in the same way as before, using f)5.16|) ()5.18[) to handle 
the asymptotics of the prefactor, then we arrive once more at (|5.32p . yielding a check on the 
rather extensive calculations. Likewise, (|5.7p gives rise to f|5.33|) . 
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5.2 Asymptotic and analytic properties of F T (?7; x, y) 

With the various representations of the function F T (rj;x,y) at our disposal, several salient fea- 
tures can be readily derived. First, it is remarkably easy to show from (|5.22j) that it has 
exponential decay for x — > — oo (as might be expected from the exponential divergence of the 
'potential' factors in the Hamiltonians H±(r];x) (|5.12p ). To be specific, we have a bound 

F T (rj; x, y) = 0(exp(aax/A)), x — > — oo. (5.35) 

Inspecting (|5.22j) . it is clear that we need only show that the integral yields a function that 
is O(l) for x —> — oo. To this end we point out that from ([A.13P we have estimates 

G(v - ia/2) = 0(exp(=Fcrat;/4), v -»• ±oo, (5.36) 

and that no poles arise for real v. Hence the function v \— > G(v — ia/2) is bounded on R. If we 
now take z — > z + (x — r\)/2 in (|5.22p . then we obtain a factor G{z — ia/2) times a factor that 
is bounded for x, y,z,r/ £ R. Thus we can invoke the bound (|5.36j) on the first factor to deduce 
that the integral is in fact bounded for x, y, r/ varying over R, completing the proof of f|5.35j) . 

It is also not hard to obtain the x — > oo asymptotics. To this end we start from the represen- 
tation (|5.32p and follow the reasoning below (|3.52p . Thus we shift the contour down by e, where 
e > is small enough so that only the simple poles at z = ±y/2 are passed (recall our standing 
assumption y > 0). The residues of the integral then follow from (|A.19p . yielding a residue sum 

/ G(±y + ia) \ 

( G R (x - rj) ) ( G ^~ V ~ ia ) ex P( ia 2/( x ~ ^)/ 2 ) + G (y ~ ia ) exp(-iay(x - rf)/2)) . (5.37) 

Using the G-asymptotics (IA.13p . it is easily verified that the remainder integral vanishes for 
x — > oo, so that we deduce 

F t (t/; x, y) ~ u T (rj; y) 1 ^ 2 exp(iaxy /2) + u T (r]] —y) 1 ^ 2 exp(—iaxy/2), x — > oo. (5.38) 
Here we have introduced the Toda u-function 

u T {rr, y) = exp(-iarjy)G(-y + ia)/G(y + ia), (5.39) 
which can also be written 

u T (.my) = c T (r/;y)/c T (r/;-y), (5.40) 
with the Toda c-function defined by 

c T (??; y) = exp(-iar]y/2)/G(y + ia). (5.41) 
The corresponding weight function is given by 

w T (y) = l/c T ( m ±y) = G(±y + ia) = As + (y)s.(y), (5.42) 

where we used the G-AAEs (|A.2p in the last step. 

The dual counterparts of these formulae are not obvious. To begin with, we have been unable 
to establish the large-y asymptotics of F T (r]; x, y). We conjecture, however, that this is given by 

F T ( ?? ; x, y) ~ e-^ 2 -^ 8 G R (x - rj) 1 / 2 exp(ia[y 2 /A + (x - V )y/2]) 

+ e lx/2+i7T/8 G R (x-ri)- 1 / 2 exp(-ia[y 2 /A + (x-ri)y/2}), y -»• oo. (?) (5.43) 
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Even when this can be shown, it is not clear whether the function Gr(x — rj) can be viewed as 
an 5-matrix for the dual dynamics. Indeed, the dual scattering theory seems quite unusual, just 
as at the classical level [17j . Moreover, like the w- function w(a — it] — iA; x + A), the u- function 
u{a — ir] — iA; x + A) has no limit for A — > oo. 
On the other hand, the similarity transforms 

Aj(r,; x) ee Gr(x - V r l/2 Hj( V ; x)G R (x - v) 1 / 2 

= exp(-ia_sd x ) + [1 + e§(— 2x - ia_s + 2r/)] ex.p(ia_sd x ), (5.44) 

can also be obtained as the limits 

Aj(rj;x) = lim Ag(a - ir] - iA;x + A), (5.45) 

A— s-oo 

cf. (|2.45p and (15. lip . Note that they have holomorphic coefficients, whereas the dual AAOs 

AJ(rj;y) = w 7 \y)- 1 ' 2 H 7 \rj;y)w 7 \y)^ 2 = ^^(exp(ia. 5 d y ) - exp(-ia_ s d y )) , (5.46) 

(which are the counterparts of A§(b; y)), have meromorphic coefficients. Surprisingly, the AAOs 
A±(a — ii] — iA; y) have no limit, whereas they do have obvious Toda counterparts, namely 

•Ajivw) = c T {ri^y 1 A 7 (ri;y)c T (ri;y) = exp{-ia_ s d y ) + ^ 1 ^ exp(m_g^) ^\ • (5.47) 

Even though w(a — %r\ — iA; x + A) has no limit either, there exists a function 

w t (t ] ;x) = 1/E(±(x-7 1 )), (5.48) 

that may be viewed as a weight function. Here, E(x) is the function featuring in (|A.20p . which 
we already had occasion to use in Section [3l cf. the paragraph containing (I3.33p . 

To explain this interpretation, we first invoke (|A.22p : This representation makes clear that 
wt{ti; x) is a real-analytic positive function on M. Secondly, we note that when we set 

ct{t}; x) = E(x - ?]), (5.49) 

then we get 

u T (m x) = c T (m x)/c T {r]; -x) = G(x - rf). (5.50) 

Hence, ignoring phases and quadratic exponentials, this it-function encodes the conjectured 
large-?/ asymptotics (|5,43p . Finally, the similarity-transformed AAOs 

Aj(x -rj) = w T (v, x)- 1/2 Hj(r]; x)w T (v, x) l/2 , (5.51) 

have holomorphic coefficients. Specifically, from (|A.23P we compute 

Aj(x) _ es(-x - ias/2) exp(—i(x + ia-$/2)K-s) 



2vr r(-i(x + ia_ 5 /2)/a s + 1/2) 



exp(ia-sd x ) + (i — > —i). (5.52) 

Thus far, we have kept x real and y positive in the function F T (ry; x, y). We proceed to study 
its analyticity features. To this end, consider the function 

U{x -rj,y) = y/a+a-WTiv, x)~ 1/2 w 7 (y)~ 1/2 F 7 '(77; x, y). (5.53) 
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The weight functions w T and wt (given by (|5,42p and (|5.48|) ) are well understood from an 
analytic viewpoint, so we need only clarify the character of H{x,y). The first point to note is 
that for each of the four integral representations (|5.22p . (|5.24p . (|5.32p and (|5.33p of F T (ri;x,y) 
the weight function factors on the right-hand side of (|5.53p ensure that the prefactors of the 
z-integrals become entire functions of x and y. Indeed, this is clear from the corresponding 
representations 

H (x, y) = E(—x) exp(3ix/2) exp ^— (y 2 — iax — a 2 /2) 

x / dzG(±z + x/2-ia/2)exp{ia[zy + z 2 /2}), (5.54) 
Jr 

H(x, y) = E{x) exp(— 3ix/2) exp {^~~^~ {y 2 + ^ ax ~ ° 2 /^) 

x / dzG(±z -x/2 -ia/2) exp {-ia[zy + z 2 /2}), (5.55) 
Jr 

H(x, y) = E(—x) exp(— ix/2) exp(ia(y 2 — x 2 )/8) 

x / dzG(z±y/2-ia)exp(-ia[zx + z 2 /2]), (5.56) 

JR+iO 

%(x, y) = E(x) exp(ix/2) exp (ia{x 2 — y 2 )/8) 

x / dzG{z ±y/2 - ia) exp (-ia[zx - z 2 /2]). (5.57) 

JR+iO 

We have singled out the function H(x,y), because it extends from the real x-axis and the 
positive y-axis (where it takes real values) to a holomorphic (i.e., entire) function in x and y. 
Taking this assertion for granted, the analytic character of F T (r]; x, y) can be read off from (|5.53p . 
We continue by proving the holomorphy claim. 

Consider first the function defined by the integral in (|5.54p . The integrand is a meromorphic 
function I(z), whose asymptotics for Rez — > ±oo readily follows from the G-asymptotics (|A.13p . 
Specifically, we get 

I(z) = 0(exp(-aRez[-Imx/2 + a/2 + Imz + Imy]), Rez— ^ oo, (5.58) 
I(z) = 0(exp(aRe z[— Imx/2 + a/2 — Imz — Imy]), Rez— >— oo. (5.59) 

Therefore, exponential decay for Re z — > oo can be achieved by taking 

Imz > Imx/2 — a/2 — Imy, (5.60) 
and for Rez — > — oo by taking 

Imz < —Imx/2 + a/2 — lmy. (5.61) 

Since the two G-functions do not depend on y, this already implies that 7i(x,y) extends to 
a holomorphic function of y. Indeed, when we continue y off the positive axis, we need only 
move the contour M up on the right and down on the left (whenever need be) so as to retain 
exponential decay. 

For the x-continuation there is also no problem coming from the tail ends of the contour, 
but we need to avoid that the contour gets pinched between the upward and downward pole 
sequences (cf. ()A.16p - ()A.17j) ). 



z = — x/2 — ia/2 — Zki, z = x/2 + ia/2 + z^, k, I £ N, 



(5.62) 
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as x is continued off the real axis. This can be achieved by requiring 

x $l —i[a, oo). (5.63) 

As a result, H(x, y) extends to a holomorphic function of x and y outside the half line (|5.63p . 

Turning to the representation (|5.55p . we can argue in the same way to conclude that expo- 
nential decay for Re z — > oo can be achieved by taking 

Imz < Imx/2 + a/2 — Imy, (5.64) 
and for Re z — > — oo by taking 

Imz > —Imx/2 — a/2 — Imy. (5.65) 
Here we get poles for 

z = x/2 — ia/2 — Zki, z = — x/2 + ia/2 + Zki, 6 N, (5.66) 

as x is continued off the real axis. Thus we should require 

x $l i[a, oo), (5.67) 

so as to avoid contour pinching. It therefore follows that H(x,y) extends to a holomorphic 
function outside the half line (|5.67[) . 

Combining these two conclusions, we deduce that H(x,y) extends to a holomorphic function 
on C 2 , as asserted. It also follows that the contour integral in (|5.54p extends to a meromorphic 
function of x and y, with poles only at x = —ia — zu- Likewise, the contour integral in (|5.55|) 
yields a meromorphic function with poles only at x = ia + zm ■ 

To conclude this account of analyticity features, we point out that the holomorphy oiH{x, y) 
can also be derived in a somewhat different way from the two representations (I5.56j) - (|5.57p . 
For these cases the two downward pole sequences at z = =Fy/2 — Zki can always be avoided by 
moving the contour up. Now, however, another type of restriction arises from the requirement of 
exponential decay on the contour tails. For the integrand in (|5.56p we need Im z < a/2 — Imx/2 
on the right tail, but to obtain exponential decay on the left tail we must require Imx > —a. 
Thus we can only deduce holomorphy for Imx > —a. Likewise, for (|5.57p we need Imx < a 
and Imz < a/2 + Imx/2 on the left tail, so we can only infer holomorphy for Imx < a. 

Even so, from these two findings we can again conclude holomorphy on C 2 . Moreover, it 
follows that the contour integrals in (|5.56p and (|5.57p give rise to meromorphic functions of x 
and y with poles only at x = — ia — Zki and x = ia + Zki, respectively. 



5.3 Joint eigenfunction properties 

To complete this section, we verify that the joint eigenfunction properties have survived the 
A — > oo limit. Due to the simpler analyticity properties of the pertinent contour integrals (com- 
pared to the hyperbolic case), this is rather straightforward. First, to show that F T (rj; x, y) is an 
eigenfunction of the Hamiltonians H±(r)\x) (|5,12p with eigenvalues 2c±(y), we need only show 
that the AAOs A±(0; x) (|5.44p have the latter eigenvalues on the function Gr(x) _1 / 2 F t '(0; x, y). 

To this end we invoke the representation (|5.26p . It follows from our analysis of the contour 
integral in (|5.55p that the contour integral in (|5.26p with rj = defines a function M(x,y) that 
extends to a meromorphic function of x and y with poles occurring solely for x = ia + zti- We 
may write this function for x ^ i[a, oo) as 



A4(x, y) = / (1zKt{x : z) exp(iazy), 
Jc 



(5.68) 
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where we have introduced 

K T (x,z) = G L (±z-x/2-ia/2). (5.69) 

Also, the choice of Imz on the horizontal tails of the contour C depends on x and y via 
(|5.64p / (|5.65|) on the right /left tail, while the tails are connected by a curve separating the 
upward and downward pole sequences (|5.66p . 
Now our task is to prove 

Aj(0; x)M(x, y) = 2c 5 (y)M(x, y). (5.70) 

By virtue of the already known meromorphy properties, it suffices to show (|5.70p for x varying 
over a rectangle Rex £ (—a, a), Imx £ (—4a, —3a) (say), while keeping y positive. The restric- 
tion on Im x ensures that the upward and downward z-pole sequences (|5.66p are at a distance at 
least 2a from the real axis. Accordingly, we choose the contour C to coincide with the real axis 
for Rez G [—2a, 2a] (say). Furthermore, we need to push down the right tail end and push up 
the left tail end sufficiently far to retain exponential decay when we act with the x-shifts under 
the integral sign. 

The crux is now that we have a kernel identity 

K T (x - ias, z) + [1 + es(-2x - ias)]K T (x + ia_$, z) 

= K T (x, z - ias/2) + K T (x, z + ias/2). (5.71) 

(This identity can be readily checked by dividing first by Kt{x, z — ias/2) and then using the 
G L -AAEs flX30p .) Thus we obtain 

Aj(0;x) [ dzK T (x,z)e iazy = [ dz(K T (x, z - ia^/2) + K T (x, z + ia_ a /2)) e iazy . (5.72) 
Jc Jc 

Shifting C up and down by a_$/2, no poles are met, and on the resulting contours C+ and C_ 
there is still exponential decay at the left and right. Hence we get 

e s (-y) [ dzK T (x,zy az y + e 5 (y)[ dzK T (x, z)e mzy . (5.73) 
Jc + Jc- 

Since the two integrands are now equal and the integrals yield the same value Ai(x,y), the 
eigenvalue property (|5.70p follows. 

We point out that the kernel identity ()5.7ip plays a role similar to the kernel identity (I3.14p 
of the hyperbolic case. In the Toda case, however, we can ensure that the z-poles stay at an 
arbitrary distance from the real axis by choosing Imx appropriately, by contrast to the v -poles 
in (|3.15p . cf. (|3.16p . Moreover, in the Toda case the decay properties on the horizontal tails of 
the contour depend on Im z, whereas the choice of Im z is irrelevant in the hyperbolic case (since 
the 'z 2 -terms' drop out in the \Kez\ — > oo asymptotics) . 

Next, we note that to prove that F T (r];x,y) is a joint eigenfunction of the dual Hamilto- 
nians H±(i]; y) (|5. 14|) with eigenvalues e-t(x), we need only show that the AAOs ^4±(0; y) (|5.46p 
have these eigenvalues on the function w T (y)~ 1 ^ 2 F T (0; x, y). 

To this end we use the representation (I5.32H . Accordingly we introduce 



M{x,y) = I dzKxix, z) exp(-iazx), (5-74) 
where 



c 



K T (y, z) = G(z ± y/2 - ia) exp(ia(y 2 /8 - z 2 /2)). (5.75) 
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Recalling our analysis of the contour integral in (|5.56p , we see that we should require first of all 
Imx > —a in (|5.74|) , Then the integral is well defined when we choose the horizontal tails of C 
equal to R (say) on the left and having Imz < a/2 — Imx/2 on the right, while the middle part is 
above the pole sequences at z = ±y/2 — z^i- Furthermore, the function Ai(x,y) is holomorphic 
in x and y for Imx > —a and extends to a meromorphic function with poles at x = —ia — Zki- 
In view of these features, we need only prove 

AJ(0; y)M(x, y) = e s (x)M(x, y), (5.76) 

for y varying over a square Rey G (—a, a), Imy G (—a, a) (say), while keeping x real. To do so, 
we choose the middle part of the contour C equal to 2ia + (—2a, 2a), and connect this part to 
(— oo, —3a) and (3a, oo) in the obvious way. Then the y-shifts can be taken under the integral 
sign without z-poles hitting the contour. 

With these analytic preliminaries in place, the key algebraic point is once more a kernel 
identity, namely, 

(K T (y + ias, z) - K T (y - ias, z)) = K T (y, z - ias/2). (5.77) 



2s s (y) 



(To check it, one need only divide by the r.h.s. and use the G-AAEs (|A.2[) .) Thus the l.h.s. 
of (|5.76p equals 

/ cLzKt{x, z — ia-s/2) exp(— 2mzx/a + a^). (5.78) 
Jc 

Shifting the contour up by a_s/2, no poles are met, and so the joint eigenvalue equations (|5.76|) 
follow. 

We conclude this subsection with some remarks. It follows from the eigenvalue features 
just proved that the holomorphic function H{x,y) is a joint eigenfunction of the four AAOs 
A±(x) (|5.52p and A5[(0;y) (|5.46p with eigenvalues 2c±(y) and e±(x), respectively. The 
coefficients of the former AAOs are entire, whereas the coefficients of the latter are mero- 
morphic. The coefficients of the AAOs ^4^(0; x) (|5.44|) are entire as well, but their joint eigen- 
function G/j(x) _1//2 F T (0; x, y) is meromorphic in x, with poles for x = ia + Zki- 

This shows by example that the meromorphic vs. entire character of the coefficients of the 
type of commuting AAO pairs at issue in this paper is compatible both with entire and with 
meromorphic joint eigenfunctions. In this connection, it should be noted that when the ratio 
a+/a_ is irrational, then the only multipliers that do no destroy the joint eigenfunction property 
are the constants. 

Another consequence worth pointing out consists of the relations 

T~L(x, has + ia^s) = %{x, has — ia^s), V (x,k,5) ECxZx {+, — }. (5.79) 
Indeed, these follow from the eigenvalue equations 

7i(x, y - ias) - U(x, y + ia_ 5 ) = 2iss(y)es(x)7i(x, y). (5.80) 

6 The nonrelativistic Toda case 

In this section we obtain the nonrelativistic counterparts of the quantities in Section [5l along the 
lines laid out in Subsection 14.21 for the hyperbolic case. Thus we switch to the parameters (|4.18p 
and momentum variable (14.19p . whereas the Toda analog of (I4.20P is the substitution 

ri = -m(PM), g>0. (6.1) 
L 1 
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For the operators H+(rj; x) ()5. 12[) and -4,+ (77; x) f)5.44j) . these substitutions entail the expansion 
Hl,Al = 2 + /3 2 # + 0(/3 4 ), 0->O, (6.2) 

where 

H = -h 2 d 2 x + fj?g 2 exp(-/ix), (6.3) 
is the nonrelativistic Toda Hamiltonian. Moreover, we clearly get 

lim HT. = lim A T _ = exp (2mu~ 1 d x ) +(«-)■ -i) =: M. (6.4) 

/9-»-0 /3->0 

In view of (|4.25|) the eigenvalue of the Hamiltonian H becomes p 2 /4, while the eigenvalue of the 
monodromy operator Ai remains 2 cosh.(irp/ hfi) . 

Turning to the dual operators, for H+(r);y) (|5.14p . vl^(r/;y) ()5.46p and JL+(r/;y) (|5.47p we 
obtain (with p > 0) 

limtfj = figp- 1/2 (exp(ihnd p ) + (i -> -i))p~ 1/2 =: F, (6.5) 

lim = i/igp" 1 ( exp(ihfid p ) — («—)■ — i)) =: A, (6.6) 

lim .AT = exp(-ihfidp) + — exp(iHud p )— =: A, (6.7) 

and their eigenvalue becomes exp(//x/2). The dual operators IZT, and ^4.^ have no limits, 
however. The operators ^4 and H are related by a similarity transformation with the square 
root of the limit function 

lim(h/3fi)- 1 w T (27r/fi,h/3;l3p/^) = w nr {p/hfi), (6.8) 

where 

w nr (£;) = 2/c sinh(7rA;) = 2ir/T(±ik), (6.9) 
and the operators A and .4 by similarity with 

lim(nPfx) 1 / 2 c T (2n/fJ,,hP,2fx- 1 ln(/3 W ); = c^g/^p/hfi), (6.10) 

/3->0 

where 

c nr (A; fc) = (2vr)^ 1/2 exp(-ifc In A)r(ifc). (6.11) 

(The limit (|6.8p is clear from (|5.42p . whereas (|6.10p follows by using (|A.25[) .) 

We proceed to obtain the nonrelativistic limit of the joint eigenfunction F T (r];x,y) for each 
of the representations f|5.25[) . (|5.26p . (|5.32[) and (|5.33p . taking x real and y positive. 

First, we have 

G R (x - V ) ^ G R (—,h(3;x + - \n (—)+ - In ( . (6.12) 
\H fi \ngj fx \hf3JJ 

Comparing this to the limit (|A.33p . we see that it applies to the (3 — > limit of (|6.12p . with the 
parameter A equal to 2. Therefore, the limit of Gr{x — rj) equals 1, a circumstance that also 
explains why we get coinciding limits for Hj and A J in (|6.2p - (|6.4p . cf. (|5.44p . 
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Next, from (|6.8p we see that the above substitutions imply 

lim — - — — — — = sinhfirp /Hu). (6.13) 
13^0 a + a- Trh 

Moreover, ay 2 clearly vanishes for /3 — > 0, so it now follows that the limits of the prefactors of 
the four representations are all equal to the square root of the right-hand side of (16.131) . 

Turning to the integrand in ()5.25|) . the plane wave exp(iazy) becomes exp(izp/h). Further- 
more, the substitutions on the two Gj^-functions yield 

G R (—,hP;±z + --- f- + -ln — +-ln — . 6.14 

\fi 2 2/x 4 (jl \ngj li \hf3JJ 

Comparing once more to the limit (jA.33[) , we see that it now applies to the /3 — > limit of (16.141) 
with A equal to 1. Therefore, a short calculation gives the limit 

exp (— 2gh~ l exp(— lix/2) cosh(//;z)) . (6.15) 

Taking z — > t/fi in the integral, we wind up with the limit function 

¥^{g/h-LLx/2,p/hLi), (6.16) 

where 

f 00 

F^ r (A;r,A;) = 2(7r~ 1 /csinh7r/c) 1/2 / dt exp(-2Ae" r cosht) cos(tife). (6.17) 

Jo 

It is readily verified that f|5.26|) also leads to the representation (|6.17j) for the limit function. 
Proceeding with the above substitutions for (I5.32h . we see that we can only get convergence of 
the exponentials for (3 — > when we first replace the integration variable z by h(3w (say) . Hence 
we get a factor h(3 up front, and the plane wave exp(— iazx) becomes exp(— ifxwx); moreover, 
the quadratic exponential converges to 1 for j3 — > 0. It therefore remains to consider the ratio 

hPexp(2iw\n{Png))/G ( — ,h{3; -h{3w ± — + - (— + hp)]. (6.18) 

\ Li 2ll 2 V A* J J 

Scaling the G-functions by li/2tt, we can invoke (|AT25|) to deduce that (pUB]) has /3 ->■ limit 

{2-kll)- 1 exp(2iw ln(g/h))T(-iw ± ip/2ft/i). (6.19) 
Putting the pieces together, we now get the limit function (I6.16|) represented as 

F^ r (A;r,A;) = -^(vr-^sinhvrA:) 172 / dtT(i(-t ± k)/2) exp(it(ln A - r)). (6.20) 

The representation (15.331) also leads to (j6.20p . Thus we obtain two different representations 
for the limit function (|6.16p . The resulting identity 

I dtexp (-2e- v cosht) cos(tfc) = ^- / dtT(i(-t ± k)/2) exp(-itv), (6.21) 

Jo 87r JR+iO 

is known and not hard to verify. Indeed, the eigenfunction transform associated with (16.161) 
amounts to the Kontorovich-Lebedev transform, and the integrals yield distinct representations 
of the modified Bessel function K ik {2e~ v ), cf. e.g. [HI 10.32.9, 10.32.13 and 10.43(v)]. 

From (|6.17|) it is obvious that 'w) nr (fc) -1 / 2 F^ r (A; r, k) extends from the positive /c-axis to an 
entire function of k. Neither from (|6.17p nor from (|6.20p entireness in r is manifest, but this is well 
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known (it can be inferred, e.g., from ODE theory). As already mentioned in the Introduction, 
the eigenfunction property for the dual operators seems not to occur in the standard sources. It 
is most easily checked for A (|6.6p by proceeding as in the relativistic case. Here it follows from 
the kernel identity 

iAT 1 + i, t) -K nr (k-i,t)) = K m (k,t-i), (6.22) 

where 

K nr (k,t) = T(-it/2±ik/2). (6.23) 

(This identity is the counterpart of (|5.77p .) 

The large-r asymptotics of F^ r (A; r, k) is well known. It is given by 

F^ r (A;r,A;)~u nr (A;A : ) 1 / 2 e irfc + n nr (A;-A : ) 1 / 2 e -^, r -> oo, (6.24) 

where 

u nr (A; k) = c nr (A; k)/c nr (X; -k), (6.25) 

and readily verified from (I6.20p . It seems much harder to obtain the large-Zc asymptotics from 
the above representations. Assuming plane-wave behavior, a consideration of the dual operators 
leads to the expectation 

F^ r (A;r, k) ~ exp(i</> + ik\nk- ik)e l{r ~ lnX)k 

+ exp(-i(l)-iklnk + ik)e- iir - lnX)k , k -> oo, (6.26) 

where (ft £ [—it, it). Indeed, just like for its relativistic counterpart (j5.43|) . this seems the simplest 
behavior that is consistent with the eigenvalues and dual eigenvalues. To be sure, for neither 
case it is a priori clear that the asymptotic behavior must involve plane waves. 

At any rate, a result pertinent to (|6.26p can be found in the literature: An asymptotic 
expansion for Ki p (x) with p > x > occurs in [32^ Section 7.13.2, formula (19)]. The dominant 
asymptotics does give rise to (I6.26P with 4> = —tt/A, but an unsettling 0(x~ 1 ) error term is 
present. (If dependence on x is included, then one would rather expect an 0(x) error term. 
Indeed, x — > corresponds to v — > oo, a limit for which the Toda potential is exponentially 
vanishing.) 



A The hyperbolic gamma function 

The hyperbolic gamma function was introduced and studied in [33] as a so-called minimal solu- 
tion of a special first order analytic difference equation. It is basically the same as Kurokawa's 
double sine |34j . Faddeev's quantum dilogarithm j35j . and Woronowicz's quantum exponential 
function [36j. (The precise connections between these functions are spelled out in Appendix A 
of our paper |37j.) In this appendix we review features of the hyperbolic gamma function 
G(a + , a_; z) that are used in the present paper; if need be, see [33] for proofs. 
Unless specified otherwise, we choose 

a+,a_ > 0, (A.l) 

and suppress the dependence of G on a + , a_. 

To begin with, G(z) can be defined as the unique minimal solution of one of the two analytic 
difference equations 

G{z + tog/2) , , , 

— : — j— = 2c- S (z), 6 = +,-, (A.2) 

G(z — ias/2) 
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that has modulus 1 for real z and satisfies G(0) = 1 (recall f|l .9j) for the notation used here); 
remarkably, this entails that the other one is then satisfied as well. It is meromorphic in z, and 
for z in the strip 

S = {z G C | \ha(z)\ < a}, (A.3) 
no poles and zeros occur. Hence we have 

G{z) = exp(ig(z)), z e S, (A.4) 

with the function g{z) being holomorphic in S. Explicitly, g(z) has the integral representation 

( \ f°° d y ( sin2y2; z \ 

g(a + ,a-;z)= —I — — — , z G S. (A.5) 

Jo V V 2smh ( a +y)sinh(a_y) a + a_y J 

From this, the following properties of the hyperbolic gamma function are immediate: 

G(— z) = 1/G(z), (reflection equation), (A. 6) 

G(a_, a + ; z) = G(a+, a_; z), (modular invariance), (A. 7) 

G(\a + , Aa_; Xz) = G(a + , a„; z), Ae(0,oo), (scale invariance), (A. 8) 



G(a + , a_; z) = G(a + , o_; — z). (A. 9) 

We have occasion to use a few more features that are less obvious, including the duplication 
formula 

G(a + ,a„;2z) = G(a+, a_; z ± ia+/4 ± ia_/4), (A. 10) 
the closely related formula 

G(a + ,2a^;2z) = G(a+, a_; z ± ia+/4), (A. 11) 
the explicit evaluation 

G(ia+/2 - ia_/2) = (a+/a_) 1/2 , (A.12) 
and the asymptotic behavior of G{z) for Re (z) — > ±oo. The latter is given by 

G(a+,a_;z) = exp (=pi ( x + az 2 /A) ) (l + 0(exp(-r|Re (z)\))) , Rc (z) -»• ±oo, (A.13) 
where the decay rate r can be any positive number satisfying 

r < a min(a + , a_), (A. 14) 
and where 

x = ^r — + — • (A-15) 



24 \o_ a 
Defining 

z k i = ika + + ila-, k, I G N = {0, 1, 2, . . .}, (A.16) 
the hyperbolic gamma function has its poles at 

z = z w z ki = ~ ia ~ z kh k,l €N, (G-poles), (A. 17) 
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and its zeros at 

z = z£ p z^=ia + Zki, k,leN, (G-zeros). (A. 18) 

The pole at — ia is simple and has residue 

lim (z + ia)G(z) = — (a+a^) 1 / 2 . (A.19) 

z—^—ia 2lT 

In view of these features, G(z) can be written as a ratio of entire functions, 

G(a + , a_; z) = _E(a + , a_; z)/E(a+, a_; —z), (A. 20) 

where E( ) has its zeros at 

z = z+, k,leN, (E-zeros). (A. 21) 

The function E(a+, a_; z) we have occasion to employ is defined in Appendix A of I; it is 
closely related to Barnes' double gamma function [38]. We need two more of its properties. 
First, from equations (A. 41) and (A. 43) in I we have 

= exp (\ ( .'r^lf , - —(e^ + e-—))) , (A.22) 

V2J y Vsinh(a+y)smh(a_y) a + a_ '// 

where 2 belongs to the strip S (|A.3p . Second, we need the AAEs it satisfies, namely, 

E(z + iax/2) i — x , . , . . 

-) r^r( = V^exp (izK 6 )/T (»*/<*-« + 1/2), * = +,-, A.23 

— las/ 2) 

= J_ln(.^y (A.24) 

cf. equations (A.46)-(A.47) in I. 

We also state two zero step size limits of the hyperbolic gamma function, which we need for 
taking nonrelativistic limits. The first one yields the relation to the Euler gamma function: 

limG(l, k;kz + i/2)exp (izln(27r/e) - ln(2vr)/2) = l/r(iz + 1/2). (A. 25) 

For the second one we need to require that z stay away from cuts given by ±i[a+/2, oo). Then 
we have 

lim -f^ ! r ^ — ' = exp ( („ _ d) In [2c + (z) ), M,deR, A.26 

uniformly on compact subsets of the cut plane. 

For the relativistic Toda setting it is expedient to employ two slightly different hyperbolic 
gamma functions defined by 

G R (z) = exp(ig R (z)), g R {z) = g{z) + X + az 2 /A, (A.27) 
G L (z) = exp{ig L {z)), g L {z) = g{z) - x - az 2 /A. (A.28) 

These functions are the unique minimal solutions of the analytic difference equations 
G R {z + ia„ s /2) 



G R (z - ia^s/2) 



l + e s (-2z), 6 = +,-, (A.29) 
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G L {z + ia_ s 2) -i , ic\ \ jr / a ori\ 

— ; — = l + e s {2z), <) = +,-, A.30 

with asymptotic behavior 

G R {z) = l + 0(exp(-r|Re(z)|)), Re(z)^±oo. (A.31) 

Furthermore, they are related by 

Gfl(;z)G£(-z) = 1. (A.32) 

The properties of the functions G R and Gl just stated are easy to infer from the corresponding 
properties of the hyperbolic gamma function. (In Appendix A of |37| we already introduced 
functions S R and Sl that differ from G R and Gl by the shift z — > z — ia.) 
Finally, we have occasion to use the limits 

]im g R (a+, a_; z±\s(a+, a-)) = { ±^ e +^2z), A = 1, } 

a_ 4.0 L { u ) A > 1, 

where 

s(a+,a_) = ^ln— , (A.34) 

which hold uniformly for z varying over arbitrary compact subsets of C. To our knowledge, 
these limits have not been obtained before. We present their proof in the next appendix. 



B Proof of (DOS) 

Since we have 

9l(z) = -g R (-z), (B.l) 

we need only show (lA.33j) for g R . Our proof actually yields a stronger result, which may be 
useful in other contexts. To state this result, we fix 



Ae(l/V2,oo), (B.2) 
and choose 5 satisfying 

<5e(l,V2), <5A>1. (B.3) 



Then we shall show 

g R (a + ,a-;z + \s(a + ,a_)) = - , - — - e+(-2z) + (^(a^" 1 ), o_ I 0, (B.4) 



2sin(7ra_/a + ) 



where the implied constant can be chosen uniformly for z varying over compact subsets of C. 
A key ingredient of our proof is the comparison function we used in Subsection III A of |33] to 
obtain the G-asymptotics ()A.13p . Specifically, we first focus on the difference 

d(a+, a-; z) = a + a-gR(a + ,a-;z) - A 2 g R (A, A; z), z £ S, (B.5) 

for the special A-choice 
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Observing that A > a (with equality for a + = a_), we deduce that d(z) is well defined and 
holomorphic for z in the strip S ()A.3p . Recalling (|1.4p . (|A. 15|) and (|A.27|) . we see that we may 
rewrite d{z) as 

d(a+,a-;z) = — [ dyl(a+,a-;y) exp(2iyz), (B.7) 
2* Jr 

where 

1 / a + a_ A 2 \ 

/(a + ,a_;y) = - ^ sinh(a+y)sinh(o _ y) " ~ T^yJ ■ 

The ^4-choice (|B.6[) is the unique one guaranteeing that 7(y) has no pole at y = 0. Specifically, 
we easily calculate 

I{y) = c(a + ,a„)y + 0(y 3 ), y -> 0, (B.9) 

where c(a+,a_) is a polynomial in a + and a_ of degree 4. 

Since we let a_ go to 0, we may and will assume from now on 

\>5—. (B.10) 

A a + y 1 

Next, we shift the y-contour up by 

r = n5/a+. (B.ll) 

On account of (|B.10|) . this ensures that only the simple pole at y = iir/a + is passed. The residue 
at this pole is readily calculated, yielding the representation 

d(z) = - . " +a ~ . e+(-2z) + p(z), z E S, (B.12) 
2 sm(7ra_/a + J 

where p is the remainder integral 

p(z) = — exp(— 2rz) [ duliu + ir) ex.p(2iuz), r = 7r5/a + , z E S. (B.13) 
2« Jr 

We are now prepared to replace z by 

z + Xs(a + ,a-), |Imz|<a+/2, ( B -14) 
so that we get, using (jB.llj) . 

a + a - a - 

d(z + As) = ; ; — re+( — 2z) 

v ; 2sin(7ra_/a + ) +v ; 

+ — a 5 j x e + (-25z) [ dul(u + m5/a + ) exp(2iu(z + As)). (B.15) 

Next, we note 

lim I(w) = — ( a -± . (b.16) 

a-io 2w \w8wh(a + w) sinlr (a + u//v2) / 

Since z is required to satisfy |Imz| < a+/2, it follows that the tt-integrand in (|B.15P remains 
bounded by a fixed L 1 (M)-function as a_ I 0. Thus we readily deduce the bound 

d( - z + Xs ^ = a ~ N e+(-2z) + 0(al A - 1 ), |Imz|<a+/2, a_ | 0, (B.17) 

a + a_ 2sm(7ra_/a+) 

with the implied constant uniform on compact subsets of the strip |Imz| < a + /2. 
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Next, we claim that we have 
A 2 

-g R (A,A;z + As) = 0(al A_1 ), \Imz\<a+/2, a_ | 0, (B.18) 



a + a_ 



uniformly on compacts of |Imz| < a+/2. Taking this claim for granted, we see from (1B.5|) and 
(|B.17[) that (|B.4[) holds true, uniformly on compacts of the latter strip. Now the AAE (|A.29|) 
with 5 = — implies 

g R (a + ,a-;z + ia + /2 + As) - g R (a + , a_;z- ia + /2 + As) 

= -iln(l + exp(-^(z + \s)\\ (B.19) 

Hence we have (using the definition (IA.34j) of s) 

g R (z + ia+/2 + \s)-g R (z-ia + /2 + \s) = o(exp(-^\n^jJ , a_ | 0, (B.20) 

uniformly on compacts of C. From this it is routine to infer that (IB.4[) also holds for z £ C, 
with the bound uniform for z varying over arbitrary C-compacts. 

It remains to prove the claim. To this end we start from the identity 

1 f'°° te~ l 

g R (A,A;z) = -b + (7rz/A), b+(w) = dt-—, (B.21) 
7T J w smh(i) 

which follows from equations (3.41)-(3.46) in |33j . Next, we note the bound 

b + (w + R) = Rexp(-2(w + R))(l + 0(l/R)), R^oo, (B.22) 

where the implied constant can be chosen uniform on C-compacts. (One can use the elementary 
integral 

/ te~ ct dt = -^e~ cx {l + cx), c> 0, (B.23) 
Jx c z 

to verify this estimate.) As a consequence, we obtain 
A 2 

g R (A,A;z + Xs) 

~ — ^a^" 1 In ( — J exp(-27r\/2z/a + ), z G C, a_ 1 0, (B.24) 

uniformly for z in C-compacts. Thus (a stronger version of) our claim follows. This concludes 
the proof of (|B~4|) . and so (fA~33|) follows as an obvious corollary. 

C Fourier transform formulas 

We fix complex numbers fj,, v satisfying 

— a < InijLi < Im v < a, a = (a+ + a_)/2. (CI) 
Hence we have 

Im (y — (i) G (0, a + + a_), (C.2) 
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and the function 

I{ji,v;x) = G{x-v)/G(x-n), (C.3) 
is pole-free in the strip 

Im v — a < Im x < Im fi-\- a. (C4) 
Also, from the G-asymptotics ()A.13p we deduce 

I(fi, u; x) = 0(exp(=FoIm (y — fi)Ke x/2)), Rex— > ±00. (C5) 

Therefore, for real y the function 

r q( x _ v \ 
F(fj,,u;y) = / dxexp(iaxy) — -, (C6) 



G(x - ii) 

is well defined, and analytic in \i and v in the region (jC.ip , Moreover, we retain exponential 
decay of the integrand when we let y vary over the strip 

|Imy| < Im(v-n)/2, (C.7) 

so F((i, v\ y) is analytic in y in this strip and given by (|C.6p . 

Our first aim is to obtain the Fourier transform (|C.6P in explicit form. Since the func- 
tion / (IC.3P has no poles in the strip (IC.4p . we can make a contour shift 

x z+ (fi + u)/2, x,zeR, (C.8) 

to deduce 



F([i, v\ y) = ex.p(iay(fi + v)/2) I dz exp(iazy)G(±z — k), 

Jr 



(C.9) 



where we have set 

n = {y-n)/2, ImKG(0,a). (CIO) 

Hence calculating F amounts to finding the cosine transform of G(±z — k) for k in the strip 
Im k £ (0, a). The reasoning in the proof of the following proposition, however, hinges on staying 
at first with (ICl6|) . 

A special case of the following result amounts to a hyperbolic analog of a trigonometric beta 
integral in Ramanujan's lost notebook. More precisely, (]C.lip with y = (/i + v)j2 amounts to 
equation (1.8) with r < in Stokman's paper [39J. Formula (IC.lip is also obtained in Chapter 5 
of van de Bult's Ph.D. thesis [40j (cf. Theorem 5.6.8 with n = 1), as a result of specializing 
more general formulas. In slightly different guises it occurred previously in various papers (the 
earliest ones being [HI H2J 03] ) , but a complete proof cannot be found there. 

Proposition C.l. The Fourier transform (|C.6p admits the explicit evaluation 

F(n, v- y) = (a+a_) 1/2 exp(ia?/(// + v)/2)G(ia - 2n)G(±y -ia + n), (C-ll) 



where fi,u satisfy (jC.ip . y satisfies (|C.7p . and k is given by (|C.10p . 
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Proof. We first take y real and choose fj,, v such that 

Im [y — fx) £ (ai, a + + o_), a\ = max(o + , a_). (C12) 

This entails that we can shift v down by ias or [i up by ia$ without leaving the region (jC.ll) , 
Hence we may shift under the integral sign and use the G-AAEs (|A.2|) to obtain 

F(n,v-ia s ;y) = 2 f dxe iaxy c_ s {x - v + ia s /2)G{x - v)/G{x - p), y G R. (C.13) 
Jr 

Recalling that F(fi, v; y) is analytic in y when (|C.7p holds, this can be rewritten as 

F(n, v - ia 8 ;y) = e_ 6 {-v + ia s /2)F(fx, v;y - ia s /2) 

+ e_ s {v - ia s /2)F(n, v; y + ia s /2). (C.14) 

Likewise, we obtain 

+ ias, v\ y) = es(-fi - ias/2)F(n, v;y - ias/2) 

+ e_ 5 (> + ia s /2)F{n, v; y + ia s /2). (C.15) 

Next, consider the contour shift x — > x + is, where 

s G (Imv — a, Im/x + a), (C.16) 

cf. (lOil) . It implies 

F([jL,v;y) = I dx ex.p(ia(x + is)y)G{x + is — v)/G{x + is — /u) 
Jr 

= exp(—asy)F(fi — is,v — is;y), y G R. (C17) 

Since ^; y) is analytic in y in the strip (f07l) . it now follows from (fCTT2|) that we have 

F(fi-is,u-is;y) = e asy F(fi,u;y), |Imy| < aj/2. (C.18) 

Now in (|C.15p we may take [i,v ^ [i — is,v — is with s G (0, Im ^ + a). Then we are entitled 
to invoke (|C.18P to get 

F(fi — is + ias, v — 2/) = e s{— ^ + is — 3ias/2)e asy F(fi, v; y — ias/2) 

+ e-s(fJ. ~is + 3ias/2)e asy F(n, v; y + ias/2). (C.19) 

In this equation we can let s converge to as, which yields 

F(fi, v - ia s ;y) = es{2y - fi - ia s /2)F(n, v;y - ias/2) 

+ e_ s (2y + A* + ia s /2)F(fi, v; y + ias/2). (C.20) 

Comparing (1C,20P and (lC.14p . we see that the difference yields a linear relation between 
F(fi, v;y-\-ias/2) and F(fi, v;y — ias/2). After some simplification, this relation can be rewritten 

as 

F([i,v;y + ias/2) c„ s {y - ia + k) 

F(/j,, v; y - ia 5 /2) c_ 5 (y + za-Kj 

Introducing 

F r (fi, v; y) = exp(ia(n + u)y/2)G(±y - ia + k), (C.22) 
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it is easy to verify that F r also satisfies this y-AAE. Since we can choose 5 = +, — and a + /a_ ^ 
Q, it readily follows that we must have 

F(ji,u;y) = C(ji,u)F r (ji,v;y), (C.23) 

with C independent of y. 

The upshot of our reasoning thus far is that when \x and v are restricted by ()C.12p . then the 
Fourier transform is of the form (|C.23p . with F r given by (|C.22p , By analyticity in /i and u, this 
relation now extends to the whole region (jC.ip and then, by analyticity in y, to the strip ()C.7p . 
In view of (|C.9p . we therefore have shown 

/ dz exp(iazy)G(±z - re) = C(/x, v)G{±y — ia + re), (C.24) 
Jr 

Imre G (0,o), |Imy| < Imre. 
It follows from this that we have 

C(jJL,u)=D(ji-u), (C.25) 
so it remains to prove 

D(fi -v) = (a + a-) 1/2 G(ia + fi-u). (C.26) 
To this end we substitute 

F(fi, v- y) = D{{jl - v) exp(ia(n + u)y/2)G(±y - ia + re) (C.27) 
in (|C.14p . and divide the result by 

D(n - v) ex.p(ia(fi + v)y/2)G(±y - ia s /2 - ia + re). (C.28) 

Using the G-AAEs, a straightforward calculation then yields 

D^-- + ias) = 
D(n - v) 

Setting 

D r (n-u) = G(n-v + ia), (C.30) 

we see that D r also satisfies the AAE ()C.29p . Thus we deduce as before 

D(p-u) = rjD r (p-v), (C.31) 

where 77 can only depend on a + and a_. 
As a result, we have now proved 



/ 



dx exp(iaxy)G(±x — re) = rjG(ia — 2re)G(±y — ia + re), (C.32) 
Imre S (0, a), |Imy| < Imre. 
To calculate 77, we choose re equal to ia_/2. Using the G-AAEs, this yields 
f 1 1 

/ dxex.p(iaxy) — — - = r]G(ia + /2 — m_/2) — — . (C.33) 
c+(a?) c_(y) 

Now the integral is elementary, yielding a+/c_(y). Using (|A.12p we then infer rj equals (a+a^) 1 / 2 , 
hence completing the proof. ■ 
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In view of (|C.32|) . the cosine transform of G(±x — k) is proportional to G(±y — ik), with 
k = ia — k. (C.34) 
To be specific, we have 

I — I / cos(axy)G(±x — n)dx = G{k — K)G(±y — k), (C.35) 
V 7T / Jo 

IniK G (0,a), y G K. 

Note that the proportionality constant can now be checked by taking the cosine transform 
of (jCX35jh 

We proceed by deriving a corollary of the above result, which we need in Subsection 15.11 
First, we fix v in the strip \m.v G (0, a), so that G{x — v) has exponential decay for x —> ±oo, 
and choose fi real, so that G(x — jj) is a phase for real x. Consider now the integral 

/ a \i/ 2 f 2 
— / (ixexp(zaxy)G(x — i/)[exp(— ix — ia/^ /4 + ia/ux/2)G(— x + //)], (C.36) 

y G E. 

In virtue of the G-asymptotics, the phase factor in square brackets converges to exp(iax 2 /4) for 
\i — > — oo. By dominated convergence, the integral therefore converges to 

/a\i/2 



/ dx exp(iaxy)G(x — v) exp(iax 2 /4), i/gl, Im^G(0,a). (C.37) 
On the other hand, by (jC.lip the integral (|C.36P equals 



exp( _ ix + ia[y(fX + u)/2 + F /4]) - ia + iz/2). (C.38) 

G(/i + y + m - i//2) 

Using once more the G-asymptotics, this has a fi — > — oo limit that can be written as 

e - i7r /4-4i XG ( y _ ia + v j 2 ^ exp (-^[(y -ia + v/2) 2 + 4(y + v/2)(ia - v) + u 2 ]^j . (C.39) 

The resulting equality of (f037l) and (lC~39l) can be rewritten in a more illuminating way by 
using the G-cousins G R COTl) and G L (TOHl) . Indeed, let us set 

z = y + u/2, 1mue(0,a), y G K, (C.40) 

so that (|C.37p and (|C.39p can be written as 

exp(— zq^ 2 /4) C — \ f dx exp(iaxz)G(x — u) exp (ia(x — ^) 2 /4) , (C41) 

and 

exp(— iav /A)e- ln/4 -' Mx G L (z - ia) exp{-iaz(ia - v)). (C.42) 

Hence we have 

/ a \ 1/2 r 

I — ) / dxexp(ia(x + ia — v)z)G{x — u) exp (ia(x — v) /4) 

= e- iw/4 - 3ix G L (z-ia), (C.43) 
with z given by (|C.40p . Choosing now u = 2is, we obtain the following corollary. 
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Corollary C.2. Letting 

w = x + ia-2is, z = y + is, sg(0,o/2), x,y £R, (C.44) 
we /lave 

exp(iawz)Gij(w — ia)dReu> = exp(— 27r/4 — 2ix)GL{z — ia), (C.45) 

where x is given by (1A.15|) . 

Thus we have recovered the Fourier transform (A. 21) in [37] . 
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